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Abstract 

Chern-Simons modified gravity is an effective extension of general relativity that 
captures leading-order, gravitational parity violation. Such an effective theory is 
motivated by anomaly cancelation in particle physics and string theory. In this 
review, we begin by providing a pedagogical derivation of the three distinct ways 
such an extension arises: (1) in particle physics, (2) from string theory and (3) 
geometrically. We then review many exact and approximate, vacuum solutions of 
the modified theory, and discuss possible matter couplings. Following this, we review 
the myriad astrophysical, solar system, gravitational wave and cosmological probes 
that bound Chern-Simons modified gravity, including discussions of cosmic baryon 
asymmetry and inflation. The review closes with a discussion of possible future 
directions in which to test and study gravitational parity violation. 
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1 Introduction 

Over the last two decades we have experienced a wealth of observational 
data in the field of cosmology and astrophysics that has been critical in guid- 
ing physicists to our current fundamental theories. The combined CMBR [TJ, 
lensing [2], large scale structure [3] and supernovae data jl] all point to an 
early universe scenario where at matter radiation equality, the universe was 
dominated by radiation and dark matter. Soon after, the universe underwent 
a transition where it became dominated by a mysterious fluid, similar to a 
cosmologcal constant (if it is not evolving), named dark energy. Both dark 
matter and dark energy might be a truly quantum gravitational effect or sim- 
ply a modification of General Relativity (GR) at large distances. 

In the absence of a full quantum gravitational theory, how do we come 
about constructing a representative effective model? One of the most impor- 
tant unifying concepts in modern physics is the gauge principle, which in fact 
played a seminal role in the unification of the strong, weak and electromag- 
netic interactions through the requirement that the action be invariant under 
a local SU(3) x SU(2) x U(l)y gauge transformation. Given the success of 
the gauge principle, much effort has been invested in various branches of the- 
oretical and mathematical physics, unveiling deep connections between gauge 
theories and geometry. The modern route to such gauge geometrical picture is 
via the interpretation of Yang-Mills gauge fields as connections on a principal 
fibre bundle and the Riemann tensor as the curvature on the tangent bun- 
dle. Although the independent combination of the Standard Model of particle 
physics and GR accounts for all four observed forces, a complete unification 
is still lacking, in spite of their semblance as gauge theories. 

Gauge principles invariably point us in the direction of a peculiar, yet 
generic modification to GR that consists of the addition of a Pontryagin or 
"Chern-Simons" (CS) term to the action. Due to its gauge principle roots, 
such an extension connects many seemingly unrelated areas of physics, in- 
cluding gravitational physics, particle physics, String Theory and Loop Quan- 
tum Gravity. This effective theory is in contrast with other GR modifications 
that are not motivated by predictive elements of a more fundamental theo- 
retical framework. One of the goals of this report is to explore the emergence 
of CS modified gravity from these theoretical lenses and to confront model- 
independent predictions with astrophysics, cosmology and particle physics ex- 
periments. 

Many have argued that studying an effective theory derived from String 
Theory or particle physics is futile because, if such a correction to the Einstein- 
Hilbert action were truly present in Nature, it would be quantum suppressed. 
In fact, String Theory does suggest that the coupling constant in front of the 
CS correction should be suppressed at least at the electroweak scale level or 
even the Planck scale level. Indeed if this were the case, the CS correction 
would be completely undetectable by any future experiments or observations. 

Other quantities exist, however, that String Theory predicts should also 
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be related to the Planck scale, yet several independent measurements and ob- 
servations suggest this is not the case. One example of this is the cosmo logical 
constant, which according to String Theory should be induced by supersym- 
metry breaking. If this breakage occurs at the electroweak scale, then the value 
of the cosmological constant should be approximately 10 45 eV 4 , while if it oc- 
curs at the Planck scale it should be about 10 112 eV 4 . We know today that 
the value of the cosmological constant is close to 10~ 3 eV 4 , no-where near the 
String Theory prediction. 

A healthy, interdisciplinary rapport has since developed between the cos- 
mology and particle communities. On the one hand, the astrophysicists con- 
tinue to make more precise and independent measurements of the cosmological 
constant. On the other hand, the particle and String Theory communities are 
now searching for new and exciting ways to explain such an observable value, 
thus pushing their models in interesting directions. 

Similarly, a healthy attitude, perhaps, is to view the CS correction as a 
model-independent avenue to investigate parity violation, its signatures and 
potential detectability, regardless of whether some models expect this correc- 
tion to be Planck suppressed. In fact, there are other models that suggest the 
CS correction could be enhanced, due to non-perturbative instanton correc- 
tions [5], interactions with fermions [6], large intrinsic curvatures [7] or small 
string couplings at late times [8 9 ,10111 1II12|I13II14II15|I16|I17| . 




After mathematically defining the effective theory [Section|2j, we shall dis- 
cuss its emergence in the Standard Model, String Theory and Loop Quantum 
Gravity [Section [3] . We begin by reviewing how the CS gravitational term 
arises from the computation of the chiral anomaly in the Standard Model 
coupled to GR. While this anomaly is cancelled in the Standard Model gauge 
group, we shall see that it persists in generic Yang-Mills gravitational theo- 
ries. We then present the Green-Schwarz anomaly canceling mechanism and 
show how CS theory arises from String Theory, leading to a Pontryagin cor- 
rection to the action in four- dimensional gravity coupled to Yang-Mills theory. 
The emergence of CS gravity in Loop Quantum Gravity is then reviewed, as 
a consequence of the scalarization of the Barbero-Immirzi parameter in the 
presence of fermions. 

Once the effective model has been introduced and motivated, we shall 
concentrate on exact and approximate solutions of the theory both in vac- 
uum and in the presence of fermions [Sections @J and E] . We shall see that 
although spherically symmetric spacetimes that are solutions in GR remain 
solutions in CS modified gravity, axially-symmetric solutions do not have the 
same fate. In the far field, we shall see that the gravitational field of a spinning 
source is CS modified, leading to a correction to frame-dragging. Moreover, the 
propagation of gravitational waves is also CS corrected, leading to an expo- 
nential enhancement /suppression of left /right-polarized waves that depends 
on wave-number, distance travelled and the entire integrated history of the 
CS coupling. 

With such exact and approximate solutions, we shall review the myri- 
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ads of astrophysical and cosmological tests of CS modified gravity [Sections [7] 
and IS]. We shall discuss Solar-System tests, including anomalous gyroscopic 
precession, which has led to the first experimental constraint of the model. We 
shall then continue with a discussion of cosmological tests, including anoma- 
lous circular polarization of the CMBR. We conclude with a brief summary of 
leptogenesis in the early universe as explained by the effective theory. 

This review paper consists of a summary of fascinating results produced 
by several different authors. Overall, we follow mostly the conventions of [T81IT9] , 
which are the same as those of [2U] and |21| . unless otherwise specified. In par- 
ticular, Latin letters at the beginning of the alphabet a,b, . . . ,h correspond to 
spacetime indices, while those at the end of the alphabet i,j,...,z stand for 
spatial indices only. Sometimes i,j,...,z will instead stand for indices repre- 
senting the angular sector in a 2 + 2 decomposition of the spacetime metric, 
but when such is the case the notation will be clear by context. Covariant 
derivatives in four (three) dimensions are denoted by V a (Dj) and partial 
derivatives by d a (di). The Levi-Civita tensor is denoted by e abcd , while e abcd is 
the Levi-Civita symbol, with convention e 0123 = +1. The notation [A] stands 
for the units of A and L stands for the unit of length, while the notation 
0(A) stands for a term of order A. Our metric signature is (— , +, +, +) and 
we shall mostly employ geometric with G = c = 1, except for a few sections 
where natural units shall be more convenient h = c = 1. 



2 The ABC of Chern-Simons and its Tools 



2. 1 Formulation 



CS modified gravity is a 4-dimensional deformation of GR, postulated 
by Jackiw and Pi [22] |_j. The modified theory can be defined in terms of its 
action: 

S '■= Seh + ScS + Stf + S'mat, (1) 

where the Einstein-Hilbert term is given by 



Si 



EH 



K / o1 4l X\ 
JV 



-gR, 



the CS term is given by 



S cs = +a~ [ d 4 xJ=gtf *RR, 
4 Jv 



the scalar field term is given by 



(2) 



(3) 



Sit 



-13- / d A Xy 

H 2 Jv 



g ab (V o 0) (V 6 tf) + 2V(&) 



(4) 



1 Similar versions of this theory were previously suggested in the context of string 
theory |23|24j . and three-dimensional topological massive gravity |25|26| . 
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and an additional, unspecified matter contribution are described by 




where £ mat is some matter Lagrangian density that does not depend on In 
these equations, = 16nG, a and (3 are dimensional coupling constants, 
g is the determinant of the metric, V a is the covariant derivative associated 
with g a b, R is the Ricci scalar, and the integrals are volume ones carried out 
everywhere on the manifold V. The quantity *RR is the Pontryagin density, 
defined via 

*D D DD *D« cd ryb I p\ 

KK.— KK— Kb Kacdi W 
where the dual Riemann-tensor is given by 

*R\ cd := l -e cdef R a bef , (7) 

with e cde f the 4-dimensional Levi-Civita tensor. Formally, *R R oc R A R, but 
here the curvature tensor is assumed to be the Riemann (torsion-free) tensor. 
We shall discuss in Sec. [6] the formulation of CS modified gravity in first-order 
form. 

Unfortunately, since its inception, CS modified gravity has been studied 
with slightly different coupling constants. We have here attempted to collect 
all ambiguities in the couplings a and f3. Depending on the dimensionality 
of a and (3, the scalar field will also have different dimensions. Let us for 
example let [a] = L A , where A is any real number. If the action is to be 
dimensionless (usually a requirement when working in natural units), it then 
follows that \&] = L~~ A , which also forces [(3\ = L 2A ~ 2 . Different sections of 
this review paper will present results with slightly different choices of these 
couplings, but such choices will be made clear at the beginning of each section. 
A common choice is a = k and (3 = 0, leading to [a] = L~ 2 and [#] = L 2 , 
which was used in [22|27ll28ll29pil^ On the other hand, 

when discussing Solar System tests of CS modified gravity, another common 
choice is a = — £/3 and (3 = — 1, where £ is some length scale associated with 
t? [37], which then implies [a] = L, [&] = L~ x and [j3] = 1 

But is there a natural choice for these coupling constants? A minimal, 
practical and tempting choice is a = 1, which then implies that d is dimen- 
sionless and that [f3] = L~ 2 , which suggests (3 oc 4fJ. From a theoretical 
standpoint, the choice of coupling constant does matter because it specifies 
the dimensions of $ and could thus affect its physical interpretation. For ex- 
ample, a coupling of the form a oc k^ 1 suggest Scs is to be thought of as 
a Planckian correction, since G = £ 2 , where l v is the Planck length. On the 
other hand, if one wishes to study the CS correction on the same footing as 

2 When working in geometrized units, a dimensionless $ can still be achieved if 
/3 ~ k, but here k is dimensionless, thus pushing all dimensions into a, which now 
possess units [a] = L 2 
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the Einstein-Hilbert term, then it is more convenient to let a = k and push all 
units into By leaving the coupling constants unspecified with a and (3 free, 
we shall be able to present generic expressions for the modified field equations, 
as well as particular results present in the literature by simply specifying the 
constants chosen in each study 

The quantity $ is the so-called CS coupling field, which is not a constant, 
but a function of spacetime, thus serving as a deformation function. Formally, 
if "d = const. CS modified gravity reduces identically to GR. This is because 
the Pontryagin term [Eq. (jBj)] can be expressed as the divergence 

V a K a = ±*RR (8) 
of the Chern-Simons topological current 

where here T is the Christoffel connection. One can now integrate <Scs by parts 
to obtain 

Scs = a (0 K a ) \ av - % [ d 4 x^ (V Q tf) K\ (10) 

2 Jv 

where the first term is usually discarded since it is evaluated on the boundary 
of the manifolcG} The second term clearly depends on the covariant derivative 
of which vanishes if i? = const, and, in that case, CS modified gravity reduces 
to GR. 

For any finite, yet arbitrarily small V a $, CS modified gravity becomes 
substantially different from GR. The quantity V a i9 can be thought of as an 
embedding coordinate, because it embeds a generalization of the standard 3- 
dimensional CS theory into a 4-dimensional spacetime. In this sense, V a "i9 and 
V a Vf>$ act as deformation parameters in the phase space of all theories. One 
can then picture GR as a stable fixed point in this phase space. Away from 
this "saddle point," CS modified gravity induces corrections to the Einstein 
equations that are proportional to the steepness of the $ deformation param- 
eter. 

The equations of motion of CS modified gravity can be obtained by vari- 
ation of the action in Eq. (CQ). Exploiting the well-known relations 

5R b acd = \7 c 5T b ad - V d 5T b ac (11) 

and 

<c = \g M (Vafydc + Vcfyad - V d 5g ac ) , (12) 

one finds 



3 The implications of discarding this boundary term will be discussed in Sec. 12.51 
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5S=k f d 4 xy^g (R ab - \g a bR + -C ab - ^-T ab ] 5g 
Jv V 2 k 2k J 

+ / d 4 x\/^g 
Jv 



(13) 



where □ := g ab V ' V& is the D'Alembertian operator and T ab is the total stress- 
energy tensor, defined via 



T 



■ab 



+ 



-g V s 9ab Sg, 



ab 



(14) 



where is the Lagrangian density of the scalar field action, it. the integrand 
Eq. (TjJ divided by \J—g, such that S$ = ^L^x. Thus, the total stress- 
energy tensor can be split into external matter contributions T^ t and a scalar 
field contribution, which is explicitly given by 



P 



(V«0) (V 6 tf) - -<? o6 (V tf) (V a tf) - g ab V{$) 



(15) 



The tensor C ab that appears in Eq. ( fl3|) is a 4- dimensional generalization 
of the 3-dimensional Cotton- York tensor, which in order to distinguish it from 
the latter we shall call the C-tensor[ji|. This quantity is given by 



G 



ab 



c€ cde(a VeR b) d + Vcd * R d(ab)c 



(16) 



where 

v a := V a tf , v ab := V a V b tf = V (a V b) tf (17) 

are the velocity and covariant acceleration of respectively. The last line of 
Eq. ffT3"j) represents surface terms that arise due to repeated integrations by 
parts. Such terms play an interesting role for the thermodynamics of black 
hole solutions, which we shall review in Sec. 12.51 

The vanishing of Eq. (fl3l) leads to the equations of motion of CS modi- 
fied gravity. The equations of motion for the metric degrees of freedom (the 
modified field equations) are simply 



G a b H C ab — ——T abj 
k 2k 



where G ab = R a b — \g a bR is the Einstein tensor. The trace-reversed form of 
the modified field equations 



Rab H Cab — — (T a b ~ ~gabT 

K 2k V 2 



(19) 



4 In the original work of [22j , the C-tensor was incorrectly called "Cotton tensor" , 
but the concept of a higher-dimensional Cotton- York tensor already exists [38] and 
differs from the definition of Eq. (|16l) 
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can be derived by noting that the C-tensor is in fact symmetric and traceless, 
where T = g ab T ab is the trace of the total stress-energy tensor. Thus, it follows 
that as in GR the modified field equations must also satisfy 

fl=-J-T = 0, (20) 
2k 

where the right-hand side holds in the absence of matter. 

The vanishing of the variation of the action also leads to an extra equation 
of motion for the CS coupling field, namely 

Pn$ = f3^-^RR, (21) 

which we recognize as the Klein-Gordon equation in the presence of a potential 
and a sourcing term. One then sees that the evolution of the CS coupling is 
not only governed by its stress-energy tensor, but also by the curvature of 
spacetime. One can in fact also derive this equation from the requirement of 
energy-momentum conservation: 

V a (G ab + C ab ) = ^-V a T ab , (22) 

2k 

where the first term on the left-hand side identically vanishes by the Bianchi 
identities, while the second is proportional to the Pontryagin density via 

V a C ab = ~v b *RR. (23) 

Equation fl2Tl) is then established from Eq. fl22|) . provided external matter de- 
grees of freedom satisfy V a T^ t = 0. Alternatively, Eq. fl22l) also tells us that, 
provided the scalar field satisfies its evolution equation [Eq. (j2T]) ]. then the 
strong equivalence principle is satisfied since matter follows geodesies deter- 
mined by the conservation of its stress-energy tensor. 

At this junction, one might worry if this set of equations is well-posed as 
an initial value problem, ie. that given generic initial data, there exists a unique 
and stable solution that is continuous on the initial data (a small change in 
the initial state leads to a small change in the final state). A restricted class of 
modified CS theories have already been shown to be well-posed as a Dirichlet 
boundary value problem [35], via the construction of a Gibbons-Hawking- 
York boundary term (see Sec. 12.51) . In principle, these results can easily be 
generalized to initial or final boundaries (Cauchy hypersurfaces), by treating 
the case where the normal vector to the boundary is timelike, and to generic 
CS field since the addition of kinetic or potential terms should not modify 
the analysis of [10132] • Such conclusions thus imply that given an initial state, 
there exists a unique final state in CS modified gravity. 

The issue well-posedness of the theory as an initial value problem, how- 
ever, remains still formally open, since the above arguments cannot necessar- 
ily be used to demonstrate that the theory is stable. Such instabilities are 
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rooted in the potential appearance of third time derivatives in the equations 
of motion. As regards to these higher-order derivatives, notice first that for 
v a = const, such derivatives do not arise and the modified field equations 
remain second-order. Moreover, notice that even for generic third time 
derivatives also vanish in a linear stability analysis, as these derivatives are 
multiplied by terms at least quadratic in the metric perturbation. Non-linear 
stabilities, however, could arise upon the full non-linear evolution of the mod- 
ified field equations, a topic that is currently being investigated. 

Formally, the CS modified equations of motion presented here [Eqs. ( 1181) 
and (|2l|) ] represent a family of theories, parameterized by the couplings a 
and (3. Of this family, two classes or formulations are particularly interesting: 
the non- dynamical framework (a arbitrary, (3 = 0) and the dynamical frame- 
work (a and (3 arbitrary but non- zero). These two formulations are actually 
two distinct theories, because in the dynamical formulation the scalar field 
introduces stress-energy into the modified field equations, which in turn forces 
vacuum spacetimes to possess a certain amount of "scalar hair." On the other 
hand, such hairy spacetimes are absent from the non-dynamical formulation, 
but this one instead acquires an additional differential constraint that might 
overconstrain it. 

One of the benefits of introducing the coupling constants a and (3 is that 
we can easily specify if we are considering the dynamical (a ^ ^ (3) or the 
non-dynamical (a ^ 0, (3 = 0) formulation. In this review article, we shall 
attempt to present as many generic expressions with a and (3 unspecified 
as possible, but when summarizing existing results we will have to focus on 
one specific formulation. On average, the non-dynamical formulation has been 
investigated much more than the fully-dynamical one, which is why the pre- 
sentation might seem slightly biased toward the non-dynamical theory. One 
should remember that this bias is not because the non-dynamical theory is 
preferred, but only because it is easier to study than the fully-dynamical sce- 
nario. It is critical then to pay close attention to the beginning of each section 
in the remainder of this review article, as we shall specify whether the results 
that are being presented correspond to the dynamical or the non-dynamical 
theory by specifying the choice of a and (3 (an issue that is particularly relevant 
when discussing solutions to the modified theory in Sees. H]and[5]). 



2.2 Non- dynamical CS gravity and the Pontryagin Constraint 

The non-dynamical framework is defined by setting (3 = at the level 
of the action, such that the scalar field does not evolve dynamically, but is 
instead externally prescribed. Such was the formulation introduced by Jackiw 
and Pi |22j, with the particular choice a = k and (3 = 0, which implies 
[0] = L 2 . 

Within this non-dynamical model, there is a particular choice of t?, pro- 



11 



posed by Jackiw and Pi [22], that has been used extensively: 

1 



■& = - -> Vf, 



0,0,0 



(24) 



where \i is some mass scale, such that \p] = L~ l . We shall refer to Eq. 
as the canonical CS coupling. This choice of CS scalar is popular because for 
certain sufficiently symmetric line elements (eg. the Schwarzschild metric), the 
4-dimensional C-tensor reduces exactly to the ordinary 3-dimensional Cotton 
tensor. Moreover, with this choice, spacetime-dependent reparameterization 
of the spatial variables and time translation remain symmetries of the CS 
modified action [22J. In spite of this, there is nothing truly "canonical" about 
this choice of embedding coordinate and other interesting choices are also 
possible. 

Irrespective of the choice non-dynamical CS modified gravity, is a con- 
strained theory, in the sense that all solutions must satisfy an additional dif- 
ferential condition, sometimes referred to as the Pontryagin constraint: 

*RR = 0. (25) 

This constraint arises directly from the variation of the action in Eq. 0131) 
with (3 = 0. We shall see in Sees. H] and [5] that this constraint imposes severe 
restrictions on the dynamics of solutions of the non-dynamical theory. 

What does the Pontryagin constraint really mean physically? Some in- 
sight can be gained by reformulating this constraint in terms of its spinorial 
decomposition. Grumiller and Yunes [53] have realized that the trivial relation 



*RR = *CC, (26) 

where C is the Weyl tensor 

C ab cd := R ab cd - 25 [ ;R d \ + ^ c 5 b d] R (27) 

and *C its dual 

*C a b cd := l -e cde fC\ e}} (28) 

opens the door to powerful spinorial methods that allows one to map the 
Weyl tensor into the Weyl spinor [UJ, which in turn can be characterized by 
the Newman-Penrose (NP) scalars (^o, ^2, ^3, ^4)- Following the nota- 
tion of [12], the Pontryagin constraint translates into a reality condition on a 
quadratic invariant of the Weyl spinor, J, 

3 (J) = 3f(* *4 + 3*l -3^3) = 0. (29) 

The reality condition of Eq. (1291) directly implies that any spacetime of 
Petrov types iV and O automatically satisfies the Pontryagin constrain, 
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while spacetimes of Petrov type D, II and I could violate it. Moreover, this 
reality condition also directly implies that not only the Kerr solution but also 
gravitational perturbations thereof violate the Pontryagin constraint. This is 
because, although ^ li3 = in this perturbed spacetime, ^(^2) 7^ 7^ ^(^2) 
generically, which violates Eq. (|29|) [35]H3] . 

Another reformulation of the Pontryagin constraint can be obtained from 
the electro-magnetic decomposition of the Weyl tensor (cf. e.g. [44]). given by 

(C abcd + l -e abef C ef cd )u h u d = E ac + iB ac , (30) 

where u a is a normalized time-like vector, and E ac and B ac are the electric and 
magnetic parts of the Weyl tensor respectively. Grumiller and Yunes [35j have 
shown that in this decomposition [15], the Pontryagin constraint reduces to 

E ab B ab = 0. (31) 

Such a restriction forces certain derivatives of the Regge- Wheeler function in 
the Regge- Wheeler [46J decomposition of the metric perturbation to vanish, 
which has drastic consequences for perturbations of the Schwarzschild space- 
time, as we shall discuss in Sec. [5] 

The electromagnetic decomposition of the Pontryagin constraint leads to 
three possible scenarios: purely electric spacetimes B ab = 0; purely magnetic 
spacetimes E ab = 0; orthogonal spacetimes, where E ab is orthogonal to B ab . 
In fact, Eq. (13T1) is a perfect analogue to the well-known electrodynamics 
condition *F F oc E • B = 0. In electromagnetism, such a condition is satisfied 
in electrostatics (B ab = 0), magnetostatics (E ab = 0) and electromagnetic 
waves (E ab B ab = 0). The Pontryagin constraint can thus be rephrased as 
"the gravitational instanton density must vanish," since the quantity *F F is 
sometimes also referred to as the "instanton density." 

The severe requirements imposed by the Pontryagin constraint on the 
space of allowed solutions, together with the arbitrariness in the choice of CS 
field, make the non-dynamical formulation rather contrived. First, different 
choices of $ will lead to sufficiently different solutions, each of these with 
different observables. Without an external prescription to decide what $ is, 
one loses the predictive power of the Einstein equations and replaces it by a 
family of possible solutions. Moreover, all choices of d so far explored are rather 
unnatural or unphysical, in that they lead to a field with infinite energy, since 
the field's kinetic energy is constant. Such fields are completely incompatible 
with the dynamical framework, which implies that results arrived at in the 
non-dynamical framework cannot be directly extended into the dynamical 
scheme. Second, the Pontryagin constraint can be thought of as a selection 
rule, that eliminates certain metrics from the space of allowed solutions. Such 
a selection rule has been found to overconstrain the modified field equations, 
to the point that only the trivial zero-solution is allowed in certain cases |34j . 

Having said this, the non-dynamical framework has been useful to qual- 
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itatively understand the effect of the CS correction on gravitational parity 
violation. Only recently has there been a serious, albeit limited effort to study 
the much more difficult dynamical formulation and preliminary results seem 
to indicate that solutions found in this framework do share many similarities 
with solutions found in the non-dynamical scheme. The non-dynamical for- 
mulation should thus be viewed as a toy-model that might help us gain some 
insight into the more realistic dynamical framework. 



2.3 Dynamical CS gravity 



The dynamical framework is defined by allowing (3 and a to be arbitrary, 
but non-zero constants. In fact, (3 cannot be assumed to be close to zero (or 
much smaller than a), because then the evolution equation for becomes 
singular. This formulation was initially introduced by Smith, et. al. j37j, with 
the particular choice a = —£/3 and (3 = —1, which implies [$] = L" 1 . In 
this model, the CS scalar field is thus not externally prescribed, but it instead 
evolves driven by the spacetime curvature. The Pontryagin constraint is then 
superseded by Eq. (l2"Tj) . which does not impose a direct and hard constraint 
on the solution space of the modified theory. Instead, it couples the evolution 
of the CS field to the modified field equations. 

The dynamical formulation, however is not completely devoid of arbitrari- 
ness. Most of this is captured in the potential V{d) that appears in Eq. 
since this is a priori unknown. In the context of string theory, the CS scalar 
field is a moduli field, which before stabilization has zero potential (ie. it 
represents a flat direction in the Calabi-Yau manifold). Stabilization of the 
moduli field occurs via supersymmetry breaking at some large energy scale, 
thus inducing an (almost incalculable) potential that is relevant only at such 
a scale. Therefore, in the string theory context, it is reasonable to neglect such 
a potential when considering classical and semi-classical scenarios. 

The arbitrariness aforementioned, however, still persists through the def- 
inition of the kinetic energy contained in the scalar field. For example, there 
is no reason to disallow scalar field actions of the form 



new f) 



- [ d^X\ 
2 Jv 



-fj 



(32) 



which leads to the following stress-energy tensor 

1 



T. 



new i9 
ab 



A l + 2/3 2 (<3tf) (V„tf)(V 6 0) 



:9ab 



^{ddf + (3 2 {ddf , (33) 



where we have used the shorthand (di3) 2 := g ab (V a t?) (V&$). The Pontryagin 
constraint would then be replaced by 



#l Q? + 2/3 2 V a (V a tf) (&&) 



-a- *RR. 
4 



(34) 



Of course, the choice of Eq. (TJJ is natural in the sense that it corresponds 
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to the Klein-Gordon action, but it should actually be the more fundamental 
theory, from which CS modified gravity is derived, that prescribes the scalar 
field Hamiltonian. In the string theory context, however, the moduli field pos- 
sess a canonical kinetic Hamiltonian, suggesting that Eq. (132j) is the correct 
prescription [27] . 

Another natural choice for the potential of the CS coupling is the ex- 
tensor itself. In other words, consider the possibility of placing the C-tensor 
on the right-hand side of the modified field equations and treating it as simply 
a non-standard stress-energy contribution. Such a possibility was studied by 
Grumiller and Yunes [32] for certain background metrics, which are solutions 
in GR but not in CS modified gravity. The Kerr metric is an example of such 
a background, for which they found that the induced C-tensor stress-energy 
violates all energy conditions. No classical matter in the observable universe 
is so far known to violate all energy conditions, thus rendering this possibility 
rather unrealistic. 



2.4 Parity Violation in CS Modified Gravity 

Precisely what type of parity violation is induced by the CS correction? 
Let us first define parity violation as the purely spatial reflection of the triad 
that defines the coordinate system. The operation P [A] = X p A is then said to 
be even, parity-preserving or symmetric when X p = +1, while it is said to be 
odd, parity- violating or antisymmetric if X p = — 1. By definition, we then have 

_ e aijk_ Note 



that P 



e aijk 



= —el, where e\ is a spatial triad, and thus P 
that parity transformations are slicing-dependent, discrete operations, where 
one must specify some spacelike hypersurface on which to operate. On the 
other hand the combined parity and time-reversal operations is a spacetime 
operation that is slicing independent. 

How does the CS modification transform under parity? First, applying 
such a transformation to the action one finds that S is invariant (ie. parity 
even) if and only if $ transforms as a pseudo-scalar P \d\ = — Applying 
such a transformation to the modified Einstein equations one finds that the 
C-tensor is invariant if and only if the covariant velocity of •& transforms as a 
vector P [v a ] = +v a , or equivalently if d IS cLS £L pseudo-scalar. 

The transformation properties of the CS scalar are not entirely free in the 
dynamical formulation. Since d must satisfy the evolution equation V a i> a cx 
*RR, we see that P [v a ] = +v a , and thus $ must be a pseudo-scalar. In the 
non-dynamical framework, however, one is free to choose •& in whichever way 
desired and thus the transformation properties of the action and field equations 
cannot be a priori determined. Of course, if one is to treat the modified theory 
as descending from string theory or particle physics, then d is required to be 
a pseudo-scalar as the dynamical theory also requires. 

Statements about the parity-transformation properties of a theory do not 
restrict the parity-properties of the solutions of the theory. A clear exam- 
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pie can be derived from Maxwell's theory of electromagnetism. Propagating 
modes (electromagnetic waves) travel at the speed of light in vacuum, but in 
the presence of a dielectric medium, they become birefringent, leading to Fara- 
day rotation. Even though the Maxwell action and field equations are clearly 
parity-preserving, solutions exist where this symmetry is not respected. An- 
other example can be obtained from GR, where the theory is clearly parity 
preserving, but solutions exist (such as the Kerr metric and certain Bianchi 
models) that do violate parity. 

Such symmetry considerations can be used to infer some properties of 
background solutions (ie. representations of the vacuum state) in the dynami- 
cal formulation. First, if one is searching for parity-symmetric solutions (as in 
the case of spherically symmetric line-elements), then *RR = 0, which forces 
9 to be constant (assuming this field has finite energy). One then sees that 
parity-even line-elements will not be CS corrected. On the other hand, if one 
is considering parity-odd spacetimes (such as the Kerr metric), then the Pon- 
tryagin density will source a non-trivial CS scalar, which will in turn modify 
the Kerr metric through the field equations. Such a correction will tend to in- 
troduce even more parity-violation in the solution, as we shall discuss further 
in Sec. EJ 

Clear signals of parity violation can be obtained by studying perturba- 
tions about the background solutions. As in the case of Maxwell theory, CS 
modified gravity has the effect of promoting the vacuum to a very special 
type of medium, in which left- and right- moving gravitational waves are en- 
hanced/suppressed with propagation distance. Such an effect is sometimes 
referred to as "amplitude birefringence," and it is analogous (but distinct) 
to electromagnetic birefringence (see Sec. [5] for a more detailed discussion of 
amplitude birefringence). The modified theory then can be said to "prefer a 
chirality," since it will tend to annihilate a certain polarization mode. 

2.5 Boundary Issues 

Opposite to common belief, GR does not admit a Dirichlet boundary- 
value problem as formulated in the previous section. This is so because the 
variation of the Ricci scalar in the Einstein-Hilbert action leads to boundary 
terms that depend both on the variation of the metric and its first normal 
derivative. In order to become a well-posed, Dirichlet boundary value prob- 
lem, the Einstein-Hilbert action must be supplemented by a boundary coun- 
terterm, so-called Gibbons-Hawking- York (GHY) term. This term cancels the 
aforementioned boundary terms thus yielding a well-posed boundary value 
problem. 

The issue of non-dynamical CS modified gravity as a well-posed boundary 
value problem has been addressed by Grumiller, et. al. [39] with the conven- 
tions a = k and (3 = 0. As in GR, the CS action as presented in Eq. ([3]) does 
not lead to a well-posed boundary value problem and counterterms must be 
added. Let us then concentrate on T> C s an d, in particular, on boundary terms 
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involving normal derivatives of the variation of the metric, neglecting irrele- 
vant terms. In [39] and in this section, "irrelevant terms" are defined as those 
that are bulk terms but not total derivatives, or as those that are boundary 
terms that vanish on the boundary. 

Let us then define the induced metric on the boundary as 

h a b ■= 9ab - n a n bl (35) 

where the boundary is a hypersurface with spacelike, outward-pointing unit 
normal n a . The extrinsic curvature is then 

K ab := h c a h d b V c n d , (36) 

which is simply the Lie derivative along n a , where V a stands for the four- 
dimensional covariant derivative operator. Note that the variation of this quan- 
tity is given by 

5K ab = h c a hyv e 5g cd , (37) 

up to irrelevant terms. 

With this machinery, one then finds that the variation of the Einstein- 
Hilbert and CS actions lead to the following boundary terms [39] : 

SS EH = -2k5 ( d 3 x Vh K, (38) 

JdV 

5Scs = -2a5 [ d 3 xVh CS(K), (39) 

JdV 

up to irrelevant terms, where one defines [39] 

CS(K) := ^ k Kl DjKki- (40) 

In Eqs. ( BSD , (J3SD and fHUl) . K := K a a is the trace of the extrinsic curvature, 
i, j, k and n stand for indices tangential and normal to the hypersurface respec- 
tively, and Di is the covariant derivative along the boundary. Equation (138]) is 
in fact the GHY term, while Eq. (|39l) is analogous to this term in CS modified 
gravity. Note that Eq. (1391) depends only on the trace-less part of the extrinsic 
curvature, and thus, it can be thought of as complementary to the GHY term. 

The boundary terms introduced upon variation of the action can be can- 
celled by addition of the following counterterms: 

S b EH — 2k d 3 xVhK, 

JdV 

S bCS = 2a [ d^xVh i? CS(K) . (41) 

JdV 

Again, Eq. (|4T!) is the GHY counterterm, while Eq. (T4T]) is a new counterterm 
required in CS modified gravity in order to guarantee a well-posed boundary 
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value problem. Interestingly, we could also have performed this analysis in 
terms of the CS current, using *RR oc v a K a . Doing so [39] 



^ J d A x^f^j -&*RR + 2j d 3 xVh §CS{K) = 

~ [ d 4 x^v a K a + [ VhVCSfr), (42) 
2 Jv JdV 



2 Jv Jav 
up to irrelevant terms, where CS^j) is given by 



CS(T) := -e ni ^T l im (d^M + -T m jp Thi ) (43) 



with T l jk the Christoffel connection. 

The CS counterterm presented above, however, only holds in an adaptive 
coordinate frame, where the lapse is set to unity and the shift vanishes. In 
covariant form, Grumiller, et. al. [39] have shown that the action 

S = K j^xy/=g (r + *R R)) 

+ 2k [ d 3 xVh (k + ^n a e abcd K b e V c K de ) 
Jav \ 2k ) 

+ a f d z xVh'J{h ab ,'d), (44) 
Jav 

has a well-posed Dirichlet boundary value problem. This action is a general- 
ization of the counterterms presented above, which holds in any frame. The 
last integral is an additional term that is intrinsic to the boundary and does 
not affect the well-posedness of the boundary value problem, yet it is essential 
for a well-defined variational principle when the boundary is pushed to spatial 
infinity [37138] . 



3 The Many Faces of Chern-Simons Gravity 



3.1 Particle Physics 



The first place we encounter the CS invariant is in the gravitational 
anomaly of the Standard Model. In this chapter, we shall give a pedagogi- 
cal review and derivation of anomalies that includes the gravitational one. 

An anomaly describes a quantum mechanical violation of a classically 
conserved current. According to Noether's theorem, invariance under a clas- 
sical continuous global symmetry group G yields the conservation of a global 
current j^, with A labelling the generators of the group G: 

d a f A = 0. (45) 

An anomaly A A is a quantum correction to the divergence of j A which renders 
it non-zero, d a j aA = A A . 
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On the one hand, gauge theories with chiral fermions usually have global 
anomalies in the chiral currents, = ^al^- Such anomalies do not lead 
to inconsistencies in the theory, but they do possess physical consequences. 
Historically, precisely this type of anomaly led to the correct prediction of 
the decay rate of pions into photons, tt q — > 77, by including the anomalous 
interaction i[Q6 abcd F ab F cd . 

On the other hand, gauge anomalies are also a statement that the quan- 
tum theory is quantum mechanically inconsistent. Gauge symmetries can be 
used to eliminate negative norm states in the quantum theory, but in order 
to remain unitary, the path integral must also remain gauge invariant. Quan- 
tum effects involving gauge interactions with fermions can spoil this gauge- 
invariance and thus lead to a loss of unitarity and render the quantum for- 
mulation inconsistent. Therefore, if one is to construct a well-defined, unitary 
quantum theory and if gauge currents are anomalous, then these anomalies 
must be cancelled by counterterms. 

A common example of a global anomaly in the Standard Model is the 
violation of the U(l) axial current by a one-loop triangle diagram between 
fermion loops and the gauge field external legs. Let us then derive the anomaly 
using Fujikawa's approach in 1 + 1 dimensions |4^ToT?j , generalized to d + 1 
dimensions. Since amplitudes and currents can be generated from the path 
integral, Fujikawa realized that anomalies arise from the non-invariance of 
the fermionic measure in the path integral under an arbitrary fermionic field 
redefinition. For concreteness let us consider a massless fermion coupled to 
electromagnetism in 3 + 1 dimensions. The action and partition function for 
this theory are the following: 

S = J d A x(-^F ah F ab + # 7 a A^) (46) 

Z = J DAD^D^e iS[Aa '^ ] , (47) 

where ip is a Dirac fermion, A a is a gauge field, F a b is the electromagnetic field 
strength tensor, e is the coupling constant or charge of the Dirac fermions and 
7 a are Dirac matrices, where the overhead bar stands for complex conjugation. 
Such a toy theory is invariant under a chiral tranformation of the form 

ip -> e ia ^ = + iaj 5 ip + . . . , (48) 

where a is a real number and 75 is the chiral Dirac matrix. Such an invari- 
ance leads to the U(l) global Noether current j„ xial = ^jajB'fp- I n the path 
integration approach, current conservation is exhibited by studying the Ward 
identities, which can be derived by requiring that the path integral be invariant 
under an arbitrary phase redefinition of the Dirac fermions. 

The non invariance of the fermionic measure is precisely the main ingre- 
dient that encodes the chiral anomaly. In order to study this effect, we must 
first define the measure precisely. For this purpose, it is helpful to expand i/j 
in terms of orthonormal eigenstates of i^ a D a : 
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l a D a ct) m = A m m (49) 

and 

^P(x) =J2 a m<Pm(x), (50) 

m 

where Grassmann variable multiplying the c-number eigenfunctions 

4> m (x) and D a is the gauge covariant derivative. The measure is then defined 

as 

DipD'tp — TT da n da n . (51) 

n 

Let us now demand invariance of the partition function: 

DADijD^e iS[a '^ ] = f DA'Dij'D^'e iS ' la ^\ (52) 



where the fermions transform as follows 

ip(x) — > ip'(x) + e(x), ^>{x) — > ty{x) + e(x) (53) 

for a chiral transformation e(x) = ia{x) r )^ ) il){x). We see then that the La- 
grangian transforms as 

d i x$'iT a D a iP') = J d^x^D^ - («9 a a)^7 a 75# (54) 

Assuming that the measure is invariant, integrating by parts and varying the 
action with respect to a, we recover the Ward identity 

d a < V7 a 75^ >= 0, (55) 

which is nothing but the statement of axial current conservation. 

Naively, we might conclude that the classical global current carries over 
to the quantum one, but Fujikawa [4~9ll5"0] realized that such a reasoning as- 
sumes that the path integral measure is invariant. A more careful analysis 
then reveals that a change of variables in the measure affects the coefficient 
of the Dirac fermion eigenstate expansion via 

a 'm = ^2($mn + B mn )a n (56) 

n 

where 

Bran = 1 J d^X <$ m (X^<\) n . (57) 

Using the Grassmanian properties of a m , this transformation returns the Ja- 
cobian in the measure 

Di)'Di)' = [det(l + B)Y 2 DipDtfj, (58) 

where det(-) and Tr(-) shall stands for the determinant and trace respectively. 
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The key to obtaining the anomaly resides in computing det(l + B). Ex- 
panding the determinant to first order in a 

det(l + B) = e^i+B)] = e Ti(B) > (59) 

and hence, 

[det(l + B)Y 2 = e -2i/rf 4 ^WE„^750n(x) (gQ) 
Regulating the fermion composite operator with a cut-off \ n /M 

E 4>n( X )lB4>n( X ) 2 (^)750n(^)e^ (61) 
n n 

and using that the mode functions are eigenfunctions of ry a .D a , we can also 
write Eq. ( IBTi) as 

E fe^0n =< \tr[^ aD ^ M2 ) > (62) 

n 

In order to simplify this expression we can use the identity (i-f a D a ) 2 = -D a D a + 
(l/2)a ab F a b where a ab = (i/2)[j a , 7 b ] which leads us to evaluate 

< x|Tr[ 75 e(- z?2+ ( 1 / 2 ^ abF -)/ M2 ]|x > . (63) 

As we take the limit M — > oo we can expand in powers of the background 
electromagnetic field by writing — D 2 = —d 2 + .... We are led to the following 
expression: 

I -d 2 /M 2 i • f d 4 kE _ k 2 /M 2 %M 2 

< x e ° > M \x >= i / -e k eI m = (64) 

1 1 J 2vr 4 16tt 2 v ; 

The other terms that will follow arise from bringing down powers of the back- 
ground field. 

Terms with one power of the background field and with the trace of 75 
vanish, since Tr[ 75 <7 a ft] = 0. In the limit M — > 00, terms that are second order 
in the background field also vanish, leaving: 



Tr 



Isl^^Fah?] < x\e- 92 ' M2 \x >= --^-e ahcd F ab F cd . (65) 



2 V 2M 2 



32tt< 



Putting all the non- vanishing terms above together gives us the Jacobian pref- 
actor: 

[det(l + B)]- 2 = e *f^M*)(T^^F ab F cd) (66) 
The partition function with this change of variables becomes 

Z[A] = J Dip Dlpe* S dH&rDai>+ a (x)(d a j?+(l/Wn 2 )e<> M F ab F cd )) ^ 

When we vary the partion function with respect to a we get the famous ABJ 
anomaly [5Tf52] 

8a Ja = -^^FabFcd. (68) 
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The above derivation of the ABJ anomaly also applies for the gravita- 
tional anomaly. Similar to Eq. (1681) . if we use the Riemann curvature ten- 
sor instead of the field strength tensor, we will obtain the gravitational ABJ 
anomaly: 

D a jt = -^-2^ abCdR ^fRc/ f . (69) 

Note that the right-hand side of this equation is proportional to the Pontryagin 
density of Eq. ([6]). The gravitational ABJ anomaly can be canceled by adding 
the appropriate counter term in the action, which in turn amounts to including 
the CS modification in the Einstein-Hilbert action. 

Recently, it has been shown that the CS action is also induced by other 
standard field theoretical means. In particular, [53] showed that the CS action 
arises through Dirac fermions couplings to a gravitational field in radiative 
fermion loop corrections, while [54.55J showed that it also arises in Yang- 
Mills theories and non-linearized gravity through the proper-time method and 
functional integration. We refer the reader to [53l54f55] for more information 
on these additional mechanisms that generate the CS action. 

3.2 String Theory 

In the previous section, we derived the chiral anomaly in a 3 + 1 gauge 
field theory coupled to fermions. We saw that while gauge and global anoma- 
lies can exist, gauge anomalies need to be cancelled to have a consistent quan- 
tum theory. In what follows we will show how the CS modification to general 
relativity arises from the Green-Schwarz anomaly canceling mechanism in het- 
erotic String Theory. The key idea is that a quantum effect due to a gauge 
field that couples to the string induces a CS term in the effective low energy 
four dimensional general relativity. 

Recall that the action of a free, one-dimensional particle can be described 
as the integral of the worldline swept out over a "target spacetime" X a (r), 
parametrized by r. The infinitesimal path length swept out is 

dl = (-ds 2 ) 1 ' 2 = (-dX a rfXV) V2 (70) 
where r] a b is a 9 + ID Minkowski target space-time and the action is then 

S = ~m J dl = -m J dr{-X a X a ) 1/2 (71) 

We can easily extend the discussion of a point particle to a string by 
parametreizing the worldsheet with a target function in terms of two coordi- 
nates (cr, r). Consider then the string, world-sheet field X a (a,r) embedded in 
a D-dimensional space-time, G a b that sweeps out a 1 + 1 world-sheet denoted 
by coordinates (a, r) and world-sheet line element ds 2 = h,ABdX A dX B . The 
indices A, B here run over world-sheet coordinates. Analogous to the point 
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particle, the free string action is described by 

S st = | J ds = | j dadrV^hh AB d a X A (a, r)d b X B (a, r)G ab , (72) 

where T is the string tension. 

This string is also charged under a U(l) symmetry and couples to the 
Neveu-Schwarz two-form potential, B^ u , via 

S B = J d 2 ad a X(a, r)d b X(a, r)B ab . (73) 

It is this fundamental string field B ab that underlies the emergence of CS 
modified gravity when String Theory is compactified to 4D. In a seminal 
work, Alvarez-Gaume and Witten [56] showed that GR in even dimensions will 
suffer from a gravitational anomaly in a manner analogous to how anomalies 
are realized in the last section. The low-energy limit of superstring theories are 
10 dimensional supergravity (SUGRA) theories. As discussed in the previous 
section, a triangle loop diagram between gravitons and fermions will generate 
a gravitational anomaly; similarly hexagon loop diagrams generate anomalies 
in 10 dimensions. Remarkably, Green and Schwarz [5"Tf58] demonstated that 
the gravitational anomaly is cancelled from a quantum effect of the string 
worldsheet B field, since the string worldsheet couples to a two form, B ab . 
The stringy quantum correction shifts the gradient of the B ab field by a CS 
3-form. This all results in modifying the three-form gauge field strength tensor 
H abc in 10D supergravity. 

H abc = {dB) abc -> {dB) abc + - (n abc (A) - a'Q abc (u)) (74) 

This naive shifting of H abc conspires to cancel the String Theory anomaly. 
We refer the interested reader to Vol II of Polchinski's book [59] for a more 
detailed discussion of the Green-Schwarz anomaly canceling mechanism. 

We begin our analysis from the compactification of the heterotic string to 
its 4D, N = 1 supergravity limit. For concreteness we consider the compact- 
ification to be on six dimensional internal space (ie. a Calabi-Yau manifold). 
Similar to the Kaluza-Klein idea, when we dimensionally reduce a 10 dimen- 
sional system to four dimensions, many fields (moduli) which characterize the 
geometry emerge. These fields cause a moduli problem since their high energy 
density will overdose the universe, hence, they need to be stabilized. The dis- 
cussion of moduli stabilization is beyond the scope of this review and we point 
the reader interested in this field of research to the work of Gukov et. al |60j . 
In what follows, we will assume that all moduli except the axion are stabilized 
and will not explicitly deal with them in our analysis. 

Our starting point is the 10D Heterotic string action in Einstein frame 
[ST] and we ignore the coupling to fermionic fields since they are not relevant 
for our discussion. In this theory the relevant bosonic field content is the 10D 
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metric, g a b, a dilaton and a set of two and three form field strength tensors 
H 3 := H abc and F 2 := F ab respectively. 



S = - ^Tr(F ab F ab " 



, (75) 



where 



H 3 = dB± 



i(n 3 (A)-aUM,) (76) 



i?2 := -B a 6, and where ^(A) := Q abc (A) and fi 3 (co>) = f2 a ;, c (u;) are the gauge 
and gravitational CS three-forms respectively, which in exterior calculus form 
are given by 

Q 3 (A) =Tr(dAAA+^AAAA.A) (77) 

We now dimensionally reduce the 10D action to 4D, N = 1 supergravity cou- 
pled to a gauge sector by choosing a four-dimensional Einstein frame metric, 
9mn = 9MN e2 ■ The 10D line element splits up into a sum of four and six 
dimensional spacetime line elements 

ds 2 10 = dsl + g mn dy m dy n , (78) 

where g mn is a fixed metric if the internal 6-dimensions are normalized to have 
volume 4a' 3 . The compactified effective gravitational action becomes 



Sad " 2k 



Ka J 



_ 2d^S*d»S 



(S + S*) 2 



(79) 



where S = + i6, with if) and 9 the four dimensional dilaton and model- 
independent axion fields respectively. The dilaton emerges as the four dimen- 
sional Yang-Mills coupling constant = e^, which we can assume here to 
be fixed, while the axion derives from the spacetime and internal components 

of B ab . 

Let us now focus our attention on the axionic sector of the 4D heterotic 
string. The bosonic low energy effective action takes the form 

S4d=-J d A x^g-(X 4 + A - ±-eT* H abc A *H abc - ]e^Tr(F ab F ab ). (80) 
a J K 12 4 

When we explicitly square the kinetic term of the three-form field strength 
tensor, 

H abc A *H abc = [(dB 2 - ^(p 3 (A) - a'Q 3 (u))] 2 (81) 

we obtain the cross term *dB 2 A Q 3 , where the dual to the three-form dB 2 is 
equivalent to exterior derivative of the axion d9 = *dB 2 . After integrating by 
parts, one ends up with the sought after gravitational Pointryagin interaction 

M 

" d 4 xf(6)RA3l (82) 
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where here f(9) = 9 V M4 P i a' and where V is a volume factor measured 
in string units and determined by the dimensionality of the compactification. 
Additionally, integration by parts also unavoidingly introduces a kinetic term 
for f(9), which we did not write explicitly above. Alexander and Gates [27] 
used this construction to place a constraint on the string scale provided that 
the gravitational CS term was responsible for inflationary leptogenesis. 

3.3 Loop Quantum Gravity 

The CS correction to the action also arises in Loop Quantum Gravity 
(LQG), which is an effort toward the quantization of GR through the postu- 
late that spacetime itself is discrete [63|64|65j . In this approach, the Einstein- 
Hilbert action is first expressed in terms of certain "connection variables" 
(essentially the connection and its conjugate momenta, the triad), such that 
it resembles Yang-Mills (YM) theory [BE] and can thus be quantized via stan- 
dard methods. Currently, there are two versions of such variables: a selfd- 
ual SL(2,C), "Ashtekar" connection, which must satisfy some reality condi- 
tions [67] ; and a real S77(2), Barbero connection, constructed to avoid the 
reality conditions of the Ashtekar one [68] . Both these formalisms can be com- 
puted from the so-called Hoist action, which consists of the Einstein-Hilbert 
term plus a new piece that depends on the dual to the curvature tensor [6U] . 
but which does not affect the equations of motion in vacuum by the Bianchi 
identities 

Ashtekar and Balachandran [70] first analyzed parity (P) and charge- 
parity (CP) conservation in LQG [70], which led them essentially to CS theory 
with a constant CS parameter. Since LQG resembles YM theory, its canonical 
variables must satisfy a Gauss-law like GR constraint D a Ej = 0, where D a is 
a covariant derivative operator and Ef is the triad. This constraint generates 
internal gauge transformations in the form of triad rotations. 

Physical observables in any quantum theory must be invariant under both 
large and small, local gauge transformations. As in YM theory, the latter can 
be associated with unitary irreducible representations of the type e me , where 
n is the winding number and 9 is an angular ambiguity parameter. Wavefunc- 
tions in the quantum theory must then be invariant under the action of these 
representations, but this generically would lead to different wavefunctions on 
different ^-sectors. Instead, one can rescale the wavefunctions to eliminate this 
9 dependance, at the cost of introducing a B-field dependence on the conju- 
gate momenta, which in turn force the Hamiltonian constraint to violate P 
and CP. 

Asthekar and Balachandran [70] noted that this ambiguity can be related, 
as in YM theory, to the possibility of adding to the Einstein-Hilbert action 
the term 



which is essentially the CS correction to the action when the scalar field $ = 9 




(83) 
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is constant and pulled out of the integral. In this sense a CS-like term arises 
naturally in LQG due to the requirement that wavefunctions, and thus physical 
observables, be invariant under large gauge transformations. 

But the ^-anomaly is not quite the same as CS modified gravity. After 
all, the above analysis is more reminiscent to the chiral anomaly in particle 
physics, discussed in Sec. 13.11 Recently, however, the connection between LQG 
and CS modified gravity has been completed, along a bit of an unexpected 
path. Taveras and Yunes [71] first investigated the possibility of promoting the 
Barbero- Immirzi (BI) parameter to a scalar field. This parameter is another 
quantization ambiguity parameter that arises in LQG and determines the 
minimum eigenvalue of the discrete area and discrete volume operators [72]. At 
a classical level, the BI parameter is a multiplicative constant that controls the 
strength of the dual curvature correction in the Hoist action [69] . Taveras and 
Yunes realized that when this parameter is promoted to a field one essentially 
recovers GR gravity in the presence of an arbitrary scalar field at a classical 
level. 

Although the Hoist action is attractive from a theoretical standpoint since 
it allows a construction of LQG in either Ashtekar or Barbero form, this 
action has also been shown to lead to torsion and parity violation when one 
couples fermions to the theory [731741175] . This issue can be corrected, while 
still allowing a mapping between GR and the Barbero- Ashtekar formalism, by 
adding to the Hoist action a torsion squared term, essentially transforming the 
Hoist term to the Nieh-Yan invariant [76]. When one couples fermions to the 
Nieh-Yan modified theory, then the resulting effective theory remains torsion 
free and parity preserving [77) . 

Inspired by the work of Taveras and Yunes [71] , Mercuri [78179) and Mer- 
curi and Taveras [80] considered the possibility of promoting the BI parameter 
to a scalar field in the Nieh-Yan corrected theory. As in the Hoist case, they 
found that the BI scalar naturally induces torsion, but this time when this 
torsion is used to construct an effective action they found that one unavoid- 
ingly obtains CS modified gravity. In particular, one recovers Eq. ([3]) with 
i? = [3/(2k)] 1//2 /3, with (3 the BI scalar and a = 3/ (327r 2 )v^3/t, while the scalar 
field action becomes Eq. (jlj) with (3 = 1 and vanishing potential. 

4 Exact Vacuum Solutions 

One of the most difficult tasks in any alternative theory of gravity is that 
of finding exact solutions, without the aid of any approximation scheme. In 
the context of string theory, Campbell, et. al. [23] showed that certain line 
elements, such as Schwarzschild and FRW, lead to an exact CS three-form, 
which thus does not affect the modified field equations. In the context of CS 
modified gravity, Jackiw and Pi [22] showed explicitly that the Schwarzschild 
metric remains a solution of the non-dynamical modified theory for the canon- 
ical choice of CS scalar. Shortly after, Guarrera and Hariton [81J showed that 
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Fig. 1. Space of solutions of Einstein gravity £ and CS modified gravity CS. In this 
figure, we have set 2k = 1 for simplicity of presentation. 

the FRW and Reissner-Nordstrom line elements also satisfy the non-dynamical 
modified field equations with the same choice of scalar, verifying the results of 
Campbell, et. al. [23]. Recently, Grumiller and Yunes [35] carried out an ex- 
tensive study of exact solutions in the non-dynamical theory for arbitrary CS 
scalars, with the hope to find one that could represent a spinning black hole. 
All of these investigations concern vacuum solutions in the non-dynamical 
framework ((3 = 0), with the coupling constant choice a = k. We shall also 
choose these conventions here. 

4-1 Classification of General Solutions 

Let us begin with a broad classification of general solutions in CS mod- 
ified gravity. Grumiller and Yunes [35] classified the space of solutions, a 2- 
dimensional representation of which is shown in Fig. HJ into an Einstein space, 
£, and a CS space, CS. Elements of the former satisfy the Einstein equations, 
while the elements of the later satisfy the CS modified field equations. The 
intersection of £ with CS, 7 := £ D CS, defines the Pontryagin space, whose 
elements satisfy both the Einstein and the CS modified field equations in- 
dependently. From the above definitions we can now classify solutions in CS 
modified gravity. Elements in 7 are GR solutions, because they satisfy the 
Einstein equations and possess a vanishing C-tensor and Pontryagin density. 
Elements in CS \ 7 are non-GR solutions, because they are not Ricci-flat but 
they do satisfy the Pontryagin constraint and the CS modified field equations. 

A full analytic study of exact solutions has been possible only regarding 
spacetimes with sufficient symmetries that allow for the modify field equations 
to simplify dramatically. For such scenarios, however, the search for CS GR 
solutions have lead mostly to either Minkowski space or the Schwarzschild 
metric. This can be perhaps understood by considering the vacuum sector of 
7, where the C-tensor becomes 

C ab \ Rah=0 = v cd *R d ^ c = v cd V d ^ c = , (84) 

where C a b c d and *C are the Weyl tensor and its dual respectively [Eqs. (1271) 
and (1281) ]. Such a condition implies the Weyl tensor must be divergenceless via 
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the contracted Bianchi identities, which leads to three distinct possibilities: 

(1) The (dual) Weyl tensor vanishes. In vacuum, elements of V are also Ricci 
flat, so this possibility leads uniquely to Minkowski space. 

(2) The covariant acceleration of d vanishes. Such a restriction imposes strong 
constraints on the geometry (cf. e.g. [12]), leading either to flat space or 
to the existence of a null Killing vector. 

(3) The contraction of the covariant acceleration and the dual Weyl tensor 
vanishes. 

Elements of 7 are very special, possessing a large number of symmetries and 
Killing vectors. On the other hand, elements of CS \ T cannot possess too 
many symmetries, which explains why it has been so difficult to find them. 

4-2 Spherically Symmetric Spacetimes 

Consider first the most general, spherically symmetric spacetime, whose 
metric can be decomposed as the warped product of two 2-dimensional met- 
rics [82.83J: a Lorentzian one g^x 1 ) (a, (3, . . . = t,r) with some coordinates 
x 7 ; and the metric on the 2-sphere Q i j(x t ) (i,j, ... = 9, 0) with some coordi- 
nates x l . Such a line element can be written in the following 2 + 2 form: 

ds 2 = g a ^) dx a dx p + <£> 2 (a; 7 ) dVt\ 2 , (85) 

where dVL\ 2 is a line element of the round 2-sphere and the warped factor is 
the square of the scalar field $(x 7 ), often called "areal" radius. 

A spherically symmetric line element [eg. Eq. (1851) ] always leads to a 
vanishing Pontryagin density, *RR = 0, and to a decoupling of the modified 
field equations [55] : 

R ab = 0, C ab = 0. (86) 
For the metric in Eq. (1551) . the only non- vanishing components of these tensors 
are R a p, R%j and C ai , provided d belongs to the generic family [55|5i] 

= F(x 7 ) + $(x 7 )G(x 4 ). (87) 

In the non-dynamical framework, these results imply that spherically sym- 
metric line elements are always pushed to CP. Similar conclusions also hold for 
spherically symmetric line-elements in non-vacuum spacetimes. 

In the dynamical framework, one must also solve the evolution equation 
for the CS scalar, which here becomes a wave-like equation without a source. 
In the absence of a potential, this wave-like equation need not necessarily have 
well-defined decaying solutions that will lead to finite energy contained in the 
scalar field. If this is the case, the scalar field is forced to be a constant, which 
reduces the modified theory to GR. 

The study of spherically symmetric line elements naturally leads to the 
study of Birkhoff 's theorem in CS modified gravity. This theorem states that 
the most general, spherically symmetric solution to the vacuum Einstein equa- 
tions is the Schwarzschild line element. For spherically symmetric line elements 
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and the CS scalar of Eq. (IHTj) . the non-dynamical CS modified field equations 
decouple and the C-tensor identically vanishes, which suffices to guarantee 
that Birkhoff's theorem still holds in the non-dynamical formulation [33] . 

In spite of the clear persistence of Birkhoff's theorem in the non-dynamical 
formulation, this theorem does not in fact hold in the dynamical framework. In 
this scheme, the dynamical field equations contain a scalar-field stress-energy 
contribution that will unavoidingly lead to non-vacuum (ie. hairy) solutions. 
Due to the presence of such a dynamical scalar field, spherical symmetry need 
not lead to staticity, and in fact, spherically symmetric spacetimes will in gen- 
eral be dynamical. Static solutions are, however, still present in dynamical CS 
modified gravity provided the CS scalar is a constant [M] . 

The study of the spherically symmetric spacetimes in CS modified gravity 
leaves us with two important lessons: 

• The existence of specific solutions depends sensitively on the choice of the 
scalar field. 

• The satisfaction of the Pontryagin constraint is a necessary but not a suffi- 
cient condition for the C-tensor to vanish. 

In fact, it is simple to construct a CS scalar, such as a trigonometric function 
of spacetime, and show that for such a scalar the C-tensor does not vanish, 
thus rendering the Schwarzschild metric not a solution of CS modified gravity, 
in spite of the vanishing of the Pontryagin density. 

4-3 Why is the Kerr Metric Not a Solution in CS Modified Gravity? 

Consider the Kerr metric in Boyer-Lindquist coordinates (t,r,6,4>): 

, 2 A - a 2 sin 2 6 , 2 AaMr sin 2 6 , , , 
ds = — at — dtdcp 

(r 2 + a 2 ) 2 - a 2 A sin 2 6 9 ^, 9 E, 9 , n7 

+ ^— — J —— sin 2 Qd<p 2 + — dr 2 + Ed6 2 (88) 

Zj za 

where E = r 2 + a 2 cos 2 and A = r 2 + a 2 — 2Mr. When CS modified gravity 
was proposed, Jackiw and Pi [22J realized that the Kerr metric would not 
be a solution of the modified theory because the Pontryagin density is not 
vanishing: 

*RR = 96^^ cos 6 (r 2 - 3a 2 cos 2 Q) (V - a 2 cos 2 0) . (89) 

This statement is also true in the dynamical frameworks, because the Pon- 
tryagin density will induce a non-constant CS scalar that will lead to a non- 
vanishing C-tensor. One can show that the parity-odd quantity in Eq. (159"]) . 
is also non- vanishing for the Kerr-Newman and Kerr-NUT spacetimes [5Tf35] . 
but it is is satisfied in certain interesting physical limits, namely a — > (the 
Schwarzschild limit) and M — > (the pp-wave limit). 

Just because the Kerr metric is not a solution in CS modified gravity does 



29 



not imply that a rotating BH solution is absent in the modified theory. Mod- 
ifications of the Kerr metric that do satisfy the CS modified field equations 
can be obtained be either considering approximate solutions or by studying 
the dynamical framework. For example, in the far field limit, M/r <C 1, the 
Pontryagin constraint is satisfied to 0(M/r) 3 and approximate modified solu- 
tions can be derived. On the other hand, in the dynamical formalism Eq. fl89l) 
serves as a source term that drives the evolution of the CS scalar, which in 
turn sources corrections in the metric [84J. 

An example of the latter can be found in the quantum-inspired studies 
of Campbell [8"5j86f24j . Reuter [87J and Kaloper [88] ■ For example, in [85|87] . 
F a b *F ab is added to the Einstein-Maxwell-Klein-Gordon action, with a 
scalar (axion) field and F a b the photon strength field tensor. The equations of 
motion for the scalar field acquire a source (the expectation value of a chiral 
current), when one treats the photon quantum mechanically. Upon accounting 
for the one- loop fluctuations of the four- vector potential [89 90J, one essentially 
finds Eq. (121]) . which on a Kerr background can be solved to find [87] 



The scalar field presents a r~ 2 fall-off and a dipolar structure, identical to 
the Kalb-Ramond axion [85]. This field leads to a non-trivial C-tensor that 
corrects the Kerr line element, thus disallowing this metric as a solution of the 
modified theory. Nonetheless, the axion in Eq. (|9"U1) could (and has) been used 
to construct correction to the Kerr metric [M], as we shall discuss in Sec. |5j 




4-4 Static and Axisymmetric Spacetimes 

Let us now consider static and axisymmetric line elements in vacuum and 
in the non-dynamical framework. Both stationary and static, axisymmetric 
spacetimes possess a timelike (d t ) a and an azimuthal (c^) a Killing vector, but 
the difference is that static metrics contain no cross-terms of type dtd<j). The 
most general such metric is diffeomorphic to [20J 



where V(p, z), Q(p,z) and A(p,z) are undetermined functions of two coordi- 
nates, p and z. 

Consider first the canonical choice of CS scalar. The modified field equa- 
tions once more decouple, as in Eq. (1861) and the spherically symmetric case, 
which implies all static and axisymmetric solutions are elements of D 5 , iden- 
tically satisfying the Pontryagin constraint. Due to this, we can make choose 
A = 1 and put the metric into Weyl class pOPS] : 




(90) 




ds 2 = -Vdt 2 + V- 1 p 2 d<f) 2 + Q 2 [dp 2 + Adz 2 ) , 



(91) 



ds i = - e w dt 2 + e -2U e ^ {dp 2 + dz 2 )+p 2 d(j) 2 ^ 



(92) 
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where U(p, z) and k(p, z) replace the functions V and fl. 
The vanishing of the Ricci tensor reduces to 



AU = 0, k p = p{U 2 p - U 2 ) , k z = 2pU p U z , (93) 

where A = d 2 / dp 2 + l/pd/d p + d 2 / dz 2 is the flat space Laplacian in cylindrical 
coordinates. The function k can be solved for through a line integral once U 
is determined |42j . 

The vanishing of the C-tensor reduces to the vanishing of the contraction 
of the dual Weyl tensor and the covariant acceleration of the CS coupling field. 
For a canonical one finds [351 



T t pt * C t(ab)p + T t zt * c t(ab)z = q ^ ( Q4 ) 

which leads to a set of nonlinear PDEs for U, in addition to the Laplace 
equation of Eq. fl93l . 

The set of equations imposed by Eq. fl86|) can be solved exactly, yielding 
either flat space, the Schwarzschild solution or the following two solutions: 

ds 2 = --dt 2 + zdz 2 + z 2 (dp 2 + p 2 d<p 2 ) (95) 

z 

and 

ds * = - ^ - l) dt 2 + " ' dz 2 + z 2 (dp 2 + sinh 2 pd<p 2 ) (96) 

These solutions, however, contain undesirable, non-physical features, such as 
the existence of a naked singularity at z = [Eq. fl95l) ]. or the existence of a 
spacelike Killing vector k a = {d t ) a in the "outside" region z > 2m [Eq. fl9"6P ]. 

The above results can also be obtained by noticing that the spatial sector 
of C a b reduces identically to the 3-dimensional Cotton- York tensor for static 
and axisymmetric line elements, which by the decoupling of the field equations 
must vanish exactly [35J . This implies the metric must be spatially conformally 
flat. Lukacs and Perjes [91] have shown that such solutions to the Einstein 
equations reduce to Minkowski, Schwarzschild or the line elements of Eqs. (1951) 
and ([9"6]) . as described above. 

The above results hold for any CS scalar in the family 

■d = -d x {t,p,z)+-d 2 {p,z,<i>). (97) 

For even more general CS scalars, as would arise for instance in the dynamical 
formulation, the modified field equations do not necessarily decouple, and 
thus, new and interesting solutions could arise. Unfortunately, when this is 
the case, the system of partial differential equations becomes too difficult to 
study analytically and has not really been analyzed in detail. 
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4-5 Stationary and Axisymmetric Spacetimes 



Consider now stationary and axisymmetric spacetimes 

ds 2 = -V (dt - wd<f)) 2 + V~ 1 p 2 d<f) 2 + Q 2 [dp 2 + Adz 2 ) , (98) 

where V, Q, A and w depend on p and z, with the latter identified with angular 
velocity of rotation about the Killing axis. 

Let us first consider the canonical choice of CS scalar, for which the 
modified field equations decouple once more. The Einstein equations can be 
used to set A to unity and put the line element in Lewis-Papapetrou-Weyl 
form H21 



ds 2 = -e 2U (dt - wd<PY + e~ 2U [e 2k (dp 2 + dz 2 ) + p 2 d(j) 2 \ , (99) 

where U and k replace V and Q. 

The vanishing of the Ricci tensor reduces to a set of partial differential 
equations similar to Eq. fl9"3"l) . with a non-trivial source that depends on w 
and an additional equation for this function. The vanishing of the C-tensor, 
however, does not lead to a vanishing Cotton- York tensor, and thus, it does 
not necessarily imply spatial conformal flatness. Grumiller and Yunes [35j have 
argued that it seems unlikely that other non-trivial and physically interesting 
solutions besides the static ones could arise, because R a b = and C a b = 
is a strongly over-constrained differential system. This is an example of how 
the non-dynamical theory can lead to an overconstrained system of equations 
when searching for physically relevant solutions, but a strict proof remains 
elusive. 

The argument presented above holds for any $ in the family 

0=- + ^, (100) 
p v 

with constants p and v. When the CS coupling is not of the form of Eq. (llOOp . 
then the field equations do not decouple, the arguments presented above do 
not hold and the system become more difficult to analyze analytically. 

In the non-dynamical framework, another route to exact solutions for 
generic CS scalars is through the Pontryagin constraint. For axisymmetric 
spacetimes, *RR ^ 0, and this leads to a complicated set of partial differential 
equations for w, U and k. Solutions of *RR = have been found to be of 
Petrov type II [55], which correspond to the Van Stockum class 



ds 2 = pVLdt 2 - 2pdtd<p + — (dp 2 + dz 2 ) , (101) 

VP v J 

where Q = Q(p, z) is arbitrary and there is no dcf 2 component. 

The reduced Van-Stockum class of metrics leads to a complete decoupling 
of the modified field equations, except for the dt 2 component. The vanishing 
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of the C-tensor can be achieved if d = z), while the dt 2 component of the 
field equations can be solved to find 



where c and d are constants. Equation (11021) is non-flat, possessing a non- 
vanishing Riemann tensor and a third Killing vector td t — <pd^ + ctd^,, but 
it satisfies R a b = and C a b = independently (ie. it belongs to D 5 ). Equa- 
tion (I103P is Riemann and Ricci non-flat (R a bcd ^ 0, R a b = —C a b ^ 0), 
belonging to CS \ 7. This last solution can be interpreted as a BH solution 
in the mathematical sense only, provided Q vanishes for some p and a Killing 
horizon appears, because it allows for closed timelike curves outside the Killing 
horizon [35J. 

We see then that physically relevant, stationary and axisymmetry, ex- 
act solutions have not been found in non-dynamical CS gravity, even when 
considering generic CS scalars fields. We say relevant solutions here, because 
non-physical ones have been found, but they either contain naked singularities 
or closed time-like curves. In fact, as an example of the latter, recently [92] 
has shown that the Godel line element ds 2 = a 2 (dt 2 — dx 2 + e 2x dy 2 /2 — dz 2 + 
2e x dtdy), which is a subclass of the metrics considered in this section, satisfies 
the non-dynamical modified field equations with the CS scalar $ = F(x, y), for 
some arbitrary function F . Such a solution is of class P because the C-tensor 
automatically vanishes with such choices. All evidence currently points at the 
modified field equations being over-constrained by the Pontryagin condition 
and the field equations. Although this evidence is strong, no proof currently 
exists to guarantee that no solution can be found. Moreover, in the dynamical 
framework, the relaxation of the Pontryagin constrain suffices to allow the 
existence of Kerr-like solutions [H] . 

4-6 PP-Waves and Boosted Black Holes 

Consider line elements that represent exact gravitational wave solutions 
(pp- waves [93]): 



where H is the only free function of u, x and y. The Aichelburg-Sexl limit [91] 
of various BHs is in fact an example of such a line element. This limit consists 
of ultrarelativistically boosting the BH, while keeping its energy finite, by 
simultaneously taking the mass to zero as the boost speed approaches that of 
light [951196] . 




For such metrics, the Pontryagin constraint vanishes identically, which at 
first seems in contradiction with the results of Jackiw and Pi [22] , who showed 
that generic linear GW perturbations lead to *RR = 0{h 2 ). In fact, these two 




(102) 
(103) 



ds 2 = —2dvdu — H(u, x, y)du 2 + dx 2 + dy 2 , 



(104) 
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results are actually consistent, because the line elements considered here and 
in [22] are intrinsically different and cannot be related via diffeomorphism. 
The fact that the Pontryagin density vanishes for the metric of Eq. (I104p is 
related to the fact that the CS velocity field v a = dj) is a null Killing vector 
v a v a = for such spacetimes [35] . 

The modified field equations decouple, except for the du 2 component. 
Consider first solutions that live in 7 [Eq. (jSBJ)], which require the vanishing 
of the C-tensor, thus forcing 

= X( u )v + ${u, x, y) . (105) 
The vanishing of the Ricci tensor forces 

AH = 0, (106) 

2_ff t yy$ \xy H X y(J) yy $ y xx) • (107) 

Therefore, for any H that solves the Laplace equation in Eq. (11061) . we can 
find a d such that Eq. (11071) is also satisfied. Such a scheme allows one to lift 
any pp-wave solution of the vacuum Einstein equations to a pp-wave solution 
of non-dynamical CS modified gravity of class 7 through a choice of $ that 
satisfies Eq. f[T05|) and flTOTj) [35]. 

Consider now solutions that live in CS \ 7. All non-uu components of 
the modified field equations decouple, and thus the vanishing of the C-tensor 
is satisfied by scalar field of the form of Eq. (I105p . Choosing X(u) = for 
simplicity, we find that the uu component of the modified field equations 
reduces to a third order PDE 

(1 + # >y d iX - #, x d !V )AH = {# >xx - $, yy )H :Xy - {H xx - H, yy )# >xy . (108) 

Simplifying this scenario further by choosing $ = a(u)x+b(u)y+c(u) , Eq. (11081) 
reduces to the Poisson equation AH = f, whose source satisfies a linear first 
order PDE bf >x — af y — f = 0, with general solution [assuming b(u) ^ 0] 

f{u, x, y) = e x/b{u) <f) [a(u)x + b(u)y] (109) 

where <fi is an arbitrary function. With such a source and two supplementary 
boundary conditions, we can now solve the Poisson equation and specify the 
full modified pp-wave solution [35] . Thus, generic CS\T solutions do exist 
in non-dynamical CS modified gravity and can be found via the algorithm 
described above. 

J^.l Non-Axisymmetric Solutions and Matter 

In the non-dynamical framework, axisymmetry seems to limit the exis- 
tence of solutions for a certain class of coupling functions. However, if either 
axisymmetry, the non-dynamical behavior of $ or the vacuum content assump- 
tion is relaxed, it is possible that solution in fact do exist. 
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Consider losing axisymmetry first. The general idea here is to add new 
degrees of freedom in the metric that could compensate for the overconstrain- 
ing of the decoupling of the modified field equations. Such idea is in fact 
inspired from approximate far-field solutions found in non-dynamical modi- 
fied CS gravity, which indeed require the presence of additional, non- vanishing, 
gravitomagnetic metric components. Up to the writing of this review, the only 
attempts to find such solutions have failed |35j . due to the incredible complex- 
ity of the differential system. 

Consider next spacetimes with matter content. The Kerr BH is a "vac- 
uum" solution in GR, but it does possess a distributional energy momentum 
tensor [97]. Moreover, in string theory and cosmology [28], CS modified grav- 
ity arises from matter currents, so the inclusion of such degrees of freedom 
might in principle be important. In the dynamical scheme, for example, one 
could lift any GR solution to a solution of CS modified gravity by requiring 
that 



where T^ at stands for the stress-energy of matter degrees of freedom (such 
as the distributional one associated with the Kerr solution), while T^ b is the 
energy-momentum of the CS coupling. One would now have solve the system of 
PDEs associated with Eq. (11 111) for the background which satisfies Eq. f II 1 D . 
Such a task, however, would imply also solving the equation of motion for 
the CS scalar, thus reducing this analysis to the study of exact solutions in 
dynamical CS modified gravity, which has not yet been performed. 

5 Approximate Vacuum Solutions 

Approximate schemes have been employed to solve the CS modified field 
equations in different limits. The first attempt along this lines was that of 
Alexander and Yunes [32|7] , who performed a far- field, PPN analysis of non- 
dynamical CS modified gravity (a = k, (3 = 0) with canonical This study 
was closely followed by that of Smith, et. al. [37J who carried out a far-field 
investigation of non-dynamical solutions representing the gravitational field 
outside a homogenous, rotating sphere, taking careful account of the match- 
ing between interior and exterior solutions. Konno, et. al. [98J investigated 
the slow-rotation limit of stationary and axisymmetric line elements in non- 
dynamical CS modified gravity with non-canonical CS scalar fields. Finally, 
gravitational wave solutions of CS modified gravity have been studied by 
a number of authors, both in Minkowski spacetime and in an FRW back- 
ground [22 28,31 29 33,36]. The last two set of studies were carried out in the 
non-dynamical formalism with the choices a = k and (3 = 0. Little is known 
about GW propagation or generation in dynamical CS modified gravity, al- 
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C ab = 8nT c 
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though early efforts are being directed on that front [Mf99] . 



5. 1 Formal Post- Newtonian Solution 

The PN approximation has seen tremendous success to model full general 
relativity in the slow-motion, weakly gravitating regime (for a recent review 
see [lOOj). This approximation is used heavily to study Solar System tests 
of alternative theories of gravity in the PPN framework, as well as to de- 
scribe gravitational waves from inspiraling compact binaries, which could be 
observed in gravitational wave detectors in the near future. For these reasons, 
it is instructive to study the PN expansion of CS modified gravity, before 
submerging ourselves in other approximate solutions. 

The PN approximation is essentially a slow-motion and weak-gravity 
scheme in which the field equations of some theory are expanded and solved 
perturbatively and iteratively. As such, this scheme makes use of multiple-scale 
perturbation theory [1U1II1U2|I1U3|I104| . where the perturbation parameters are 
the self-gravity of the objects (an expansion in powers of Newton's constant 
G) and their typical velocities v (an expansion in inverse powers of the speed of 
light). For example, matter densities p are dominant over pressures p and spe- 
cific energy densities II, while spatial derivatives are dominant over temporal 
ones. 

The PN framework also requires the presence of external matter degrees 
of freedom, ie. bodies that are self-gravitating and slowly-moving. Such objects 
can be described in a point-particle approximation |100j , or alternatively with 
a perfect fluid stress-energy tensor [105J: 



where u a is the object's four- velocity. In GR, the internal structure of the 
gravitating objects can be neglected to rather large PN order |106j . and thus 
one can effectively take the radius of the fluid balls to zero, which reproduces 
the point-particle result. This statement is that of the effacing principle [106] , 
which is the view we shall take in the next section when we study CS modified 
gravity in the PPN framework. However, care must be taken, since the effacing 
principle need not hold in alternative theories of gravity. In fact, as we shall see 
later on, the effacing principle must be corrected in CS modified gravity due to 
modifications to the junction conditions [I07I1U8I1U9I11UI111I112ITI5I2IIII3] . 




Perturbation theory, and thus the PN approximation, requires the use of a 
specific background and coordinate system. In the traditional PN scheme, one 
linearizes the field equation with the metric g ab = i] ab + h ab , where p ab is the 
Minkowski background, since cosmological effects are usually subdominant. 
Moreover, a Lorentz gauge is usually chosen h ba , a = h jb /2, which allows one 
to cast the field equations as a wave equation with non-trivial, non-linear 
source terms. One can show that to first order in the metric perturbation, the 
linearized CS modified field equations in the non-dynamical formalism, with 



T 



\ab 



(p + pll + p) u a u b + pg ab , 



(112) 
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canonical CS scalar and in the Lorentz gauge, can be written as [32f7] 

U v -K ah = -16tt (T ab - ^g ab T^j + 0(h) 2 , (113) 

where the superpotential ^K ab captures the CS modification to the PN expan- 
sion to the modified field equations and it is given by 

?U := h ab + ■& e° c \ah b)d , c , (H4) 

where e abcd stands for the Levi-Civita tensor density. The 00 component, the 
symmetric spatial part and the trace of the superpotential are equal to that 
of the metric perturbation, because the Levi-Civita symbol forces the CS cor- 
rection to vanish. The formal solution to the modified field equations then 
reduces to 

■K ab = -16tt n- 1 (r ab - ± 9ab T) + o{h)\ (115) 

where the inverse D'Alembertian operator stands for a Green function integral. 
This formal solution is in fact identical to that of the PN expansion of GR in 
the limit i? — ► 0. 

This formal solution can also be cast into a more practical form by per- 
turbatively solving for the metric perturbation [32|I7] . Let us then make the 
ansatz 

h ab = hf b K) +$( ab + 0(h) 2 , (116) 
where h ab is the GR solution, which is ^-independent and satisfies 

hf b R) := -167T D- 1 (T ab - ^g ab T^j , (117) 

and where ( ab is an unknown function that is first order in When we combine 
Eqs. (Tim (HUD and (JESSD we find 

U + $i 0cd (a(b) d ,c = 167re° c V^" 1 (T b)d - ±g b)d T) , (118) 

which can be solved for to find 

U = lGKfc^n- 1 (r b]d - ^g b)d T^j , (119) 

where we have neglected terms second order in The metric perturbation 
then reduces to 

h ab = -16tt D" 1 (T ab - \riabT) + Wn^D- 1 (d i(a T b)e)k - U l{a r] b)t T^ , 

(120) 

where k stand for spatial indices only. 
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5.2 Parameterized Post- Newtonian Expansion 



Solar system tests of alternative theories of gravity are best performed 
within the PPN framework. This framework was first proposed by Eddington, 
Robertson and Schiff |115|fl l6j. but it matured with the work of Nordtvedt 
and Will |117|I118|I119II120] (for a review see e. g. |115j ). PPN theory proposes 
the construction of a model-independent super-metric that represents the PN 
approximate solution to a family of gravity theories, parameterized by PPN 
parameters. Solar system experiments can then measure these parameters, 
thus selecting a particular member of this family. Currently, many of these 
parameters have been experimentally determined with tight error bars, all of 
which are consistent with GR [1211 . 

The PPN framework allows for tests of alternative theories of gravity 
through such PPN parameters. Given an alternative theory, one must first con- 
struct its PN solution and then compare it to the PPN super-metric. Through 
this comparison, one can read off how the PPN parameters depend on funda- 
mental parameters of the alternative theory. But since PPN parameters have 
been experimental constrained, one can propagate these constraints to the 
fundamental parameters of the alternative theory under consideration, thus 
obtaining an automatic Solar System test. 

Solar System tests in the PPN framework require the PN solution to the 
modified field equations. In this framework, however, it is not sufficient to 
leave the solution expressed in terms of the inverse D'Alembertian operator, 
but instead it must be parameterized in terms of PPN potentials, which are 
simply Green function integrals over the stress energy tensor. Moreover, the 
PN expansion must be carried out to slightly different orders in v for different 
components of the metric, so as to obtain a consistent Lagrangian formulation 
of the theoryt!] • In the PPN framework, it suffices to compute goo to 0(t> 4 ), 
g i to 0(v 3 ) and g, t j to 0(v 2 ). Finally, the Lorentz gauge differs slightly from 
the PPN gauge, related via an infinitesimal gauge transformation, where the 
latter is perturbatively defined via |115j 



where k stand for spatial indices only in the remaining of this section, h k k 
is the spatial trace of the metric perturbation and the symbol 0(A) stands for 
terms of 0(e A ), with e the standard PN expansion parameter of 0(1/ c) [32T7] . 

With this machinery at hand, we can now perturbatively expand the 
trace-reversed CS modified field equations. For the remaining of this chapter, 
we shall concentrate on the non-dynamical formalism of CS modified gravity, 
with the canonical choice of CS scalar. Expressions for the linearized Ricci 
and C-tensors to 0(4) in an arbitrary gauge and for generic i? are long and 

5 For example, such order counting is necessary in order to calculate the gravita- 
tional deflection of light consistently to first order. 



h jk k -\h^0(A) 




fc,o = 0(5) 



(121) 
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unilluminating, so we shall not present them here, but they can be found 
in [37]. On the other hand, in the PPN gauge and with the canonical choice 
for the linearized Ricci tensor becomes 

Roo = -^V 2 /i 00 - \h m V* + \h ij h m>ij + 0(6), 
Roi = — 2 V 2 /ioi — j^oo.oi + 0(5), 

% = -^V 2 % + 0(4), (122) 

while the linearized C-tensor is given by 
Coo = 0(6), 

C 0l = -^i m t V 2 h 0l , k + O(5), 

Cij = -\^^%%k + 0(4), (123) 

where V = rj^didj is the flat-space Laplacian. As in Sec. 15.11 we find two 
distinct corrections due to the CS modification: one to the transverse-traceless 
part of the spatial metric and the other to the vector metric perturbation. 

Let us now solve the linearized CS modified field equations iteratively 
and perturbatively. To 0(2), the 00 component of the metric is not modified 
by the C-tensor and the field equation becomes 

V 2 /i 00 = -8tt P , (124) 

because T = —p. Equation (11241) is the Poisson equation, whose solution in 
terms of PPN potentials is 

/ioo = 2E/ + 0(4), (125) 

where U is the Newtonian potential [32f7j . We shall not present these poten- 
tials here, but they can be found in [115] . 

To this same order, the Oi component does not provide any information, 
while the ij sector leads to the following field equations: 

V 2 ^. + tfe%V\)^ = -torpdij, (126) 

where 5{j is the Kronecker delta. This equation can be rewritten in terms of 
the superpotential of Sec. 15.11 as 

V 2 % = -87Tp5 i:i , (127) 

whose solution is [3"2T?] 

% = 2U6ij + 0(4). (128) 
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Following the same procedure as in Sec. 15.11 we can use the decomposition of 
Eq. (pID to find that 

(i j + #<? k \i( j) i, k = 0. (129) 
The second term on the left-hand side of Eq. (11291) is a second order correction 
in $ and can thus be neglected, which then renders hij = 'Kij to 0(2), where 
the latter is given in Eq. (I128p . The physical reason why the CS modification 
does not correct the spatial sector of the metric is related to the source studied 
here, together with the PPN gauge. In fact, we shall see in Sec. 15.61 that when 
one studies gravitational wave propagation in vacuum, the spatial sector of the 
metric is indeed modified. Moreover, if we were to study the 0(4) corrections 
to the spatial sector of the metric, we would probably find non-vanishing CS 
corrections, but such a study has not yet been carried out. 

To next order, 0(3), the only relevant field equations are related to the 
gravitomagnetic sector of the metric. The field equations become [3"2"t?] 

V 2 /i 0i + ^ Voi + 7^e° H *V 2 Vfc = IStt/w,, (130) 

which with the lower-order solution of Eq. (I125P and the superpotential of 
Eq. (11141) becomes 

V 2 ^ 0i + U 0i = Wir/Wi, (131) 

whose solution is 

^» = -\Vi-\w i , (132) 

where Vi and Wi are PPN vector potentials. Combining Eq. (11141) with Eq. (11161) . 
we can solve Eq. fTl32|) to find plTT] 

h oi = --Vi - l -Wi + 2d (V x V), + 0(5), (133) 

where (V x A) 1 = e y djAf. is the standard curl operator of flat space 

Finally, to next order, 0(4), we need only analyze the 00 component of 
the modified field equations. Since the 00 component of the C-tensor does not 
contribute to the field equations to this order, the modified field equations 
reduce exactly to those of GR, as also does their solution. We therefore find 
that the PPN solution to the CS modified field equations is given by 

^00 = ^ + 0(6), 

^ = ^ + 2^x^ + 0(5), 

^■ = ^ + 0(4), (134) 

(C 1 Ft) 

where g ab is the PPN solution of GR. This solution in fact satisfies the 
Pontryagin constraint *RR = to leading order because the contraction of the 
Levi-Civita symbol with two partial derivatives vanishes. Thus, the equations 
of motion for the fluid can be obtained directly from the covariant derivative 
of the stress-energy tensor (the strong equivalence principle holds). 
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By comparing this PPN solution to the super-metric of the PPN frame- 
work, we can now read off the PPN parameters of CS gravity [32|7] . Doing so, 
one finds that all PPN parameters of CS gravity are identical to those of GR 
(7 = (3 = 1, ( = and «i = a 2 = a 3 = £1 = £ 2 = £3 = £4 = 0), except for the 
new term in g^. This CS correction cannot be mimicked by any standard PPN 
parameter, and thus, in order to model parity- violating theories an additional 
PPN parameter \ must be introduced, namely 

g 0t = g£ PN) + x(MVxV) i , (135) 

where we have multiplied the curl operator by the total mass M, in order 
to make \ a proper dimensionless parameter. In CS modified gravity, this 
parameter is simply \ = 2d/M. 

Canonical CS modified gravity can then be tested by experimentally con- 
straining x, which thus directly places a constraint on the canonical CS scalar 
1? = 1/ li. The only requirement for such a contribution to be non- vanishing is 
that the PPN vector potential V % be non-vanishing. This is satisfied by non- 
static sources, ie. objects that are either moving or spinning relative to the 
PPN rest frame. Just because the correction to the gravitomagnetic sector of 
the metric is non- vanishing, however, does not imply that physical observables, 
such as the Lense-Thirring effect, will also be corrected. In fact, as we shall 
see in Sec. I7.1[ for such corrections to emerge, the curl of the gravitomagnetic 
correction to the metric must be non-vanishing, which holds only for moving 
point particles. 

A caveat should be discussed at this point. As we already alluded to, 
the point-particle approximation holds in GR for many sources, such as black 
holes, neutron stars or other regular stars, because the effacing principle holds 
and Birkhoff 's theorem also holds. In CS gravity, however, the former does not 
necessarily hold and the junction conditions are modified. One should then also 
be aware that a homogenous solution to the linearized field equations might 
be lacking from this analysis. Such a solution is most likely oscillatory in 
nature and would not affect the average behavior of the correction, probably 
acquiring importance only when the latter vanishes. In the next section we 
shall see how such a boundary solution arises and contributes significantly for 
spinning bodies, for which the homogeneous correction to the Lense-Thirring 
effect vanishes. 

5.3 Rotating Extended Bodies and the Gravitomagnetic Analogy 

The analysis presented so far has dealt with the non-dynamical formal- 
ism, concentrating primarily on point particle sources. Approximate solutions, 
however, can be developed for extended objects as well. As already described, 
in GR the effacing principle guarantees that for certain sources both the ex- 
tended and point particle approach coincide up to rather high PN order. This 
result somewhat relies on the junction conditions of GR, which here are mod- 
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ified due to the presence of the C-tensor in the field equations. Moreover, the 
PN scheme is not the only formalism in which far field solution can be in- 
vestigated. In particular, the gravitoelectromagnetic analogy |122lll23j is also 
useful to obtain far field solutions, since it allows one to use machinery from 
the theory of electromagnetism. 

Smith, et. al. [37] performed the first detailed analysis of rotating, ex- 
tended objects in CS modified gravity in the gravitoelectromagnetic scheme, 
modeled via the stress-energy tensor 

Tab = ^t( a Jb) — pt a tb, (136) 

where the current J a := (— p, J) and t a = [1,0,0,0] is a time-like unit vec- 
toi0. Such a choice of stress-energy is equivalent ot that of a pressureless 
perfect fluid in the limit as the fluid balls tend to particles. Moreover, Smith, 
et. al. [37J focused on the dynamical formulation of CS modified gravity, with 
the conventions a = —£/3 and (3 = —1. Such a study is forbidden in the 
non-dynamical formalism, because the Pontryagin constraint does not allow 
for rotating solutions. In the dynamical scheme, however, such solutions are 
allowed, provided the Pontryagin density is balanced by the dynamics of the 
scalar field. 

In spite of working in the dynamical theory, Smith, et. al. [37] chose 
a canonical thus implicitly treating this field as non-dynamical. Such a 
choice of scalar field is formally inconsistent with dynamical CS modified grav- 
ity, since the non-vanishing Pontryagin constraint will force spatial variations 
on Nonetheless, this choice could become a good approximation, if one is 
only concerned with far-field solutions. This is because, as we have shown in 
Eq. (189]) of Sec. l4.3[ the exterior Pontryagin density for a rotating source scales 
as *RR^t: aM 2 /r 7 , which then forces 

/O? = -24a^£ = 8^, (137) 

with uncontrolled remainders of 0(M/r) 4 . Thus, spatial variations will induce 
modifications to the canonical •& that are at most of 0(M/r) 3 , which is beyond 
the order considered in this section. A full dynamical study of this problem 
has yet to be carried out. 

Let us now return to the gravitoelectromagnetic formalism, in which cer- 
tain components of the metric perturbation are identified with a four-vector 
potential, namely 

A a := -h ab t b , (138) 

where h a b is the trace-reversed metric perturbation h ao := h a b — ^f] a bh and 
h = t] ab h a b is the trace of the latter. In analogy with electromagnetism, one 

6 Recall that indices are raised and lowered with the Minkowski metric in the far- 
field approximation. 
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can now define electric and magnetic fields via 

E l = diA - d A, B l = e^djAk, (139) 

where i,j stand for spatial indices only. The geodesic equations can then be 
written in terms of these fields as 

F i — —mE i — 4m (v x Bf , (140) 

where the cross product is that of flat Euclidean space, from which one can 
read off the Lorentz force acceleration a 1 = F l /m. 

We can now apply this formalism to the linearized CS modified field equa- 
tions. In the Lorentz gauge, A a a = 0, the modified field equations become 



V-5 = 0, (141) 
VxE = - — , (142) 

V-E = 4nG(p + T° ) ) (143) 
dE 1 

VxB-- □ J B = 4vrGJ, (144) 

at m cs 

where m cs is a characteristic mass scale defined via 

m cs := — r = r. (145) 

oe& M K J 

In the above equations, the first two arise from the definitions of the four- 
vector potential, while the last two (Gauss's law in Eq. ( 1143jl and Ampere's 
law in Eq. ( 1144R ) are modified in CS modified gravity, the former arising from 
the time-time component of the modified field equation and the latter from the 
gravitomagnetic sector. Smith, et. al. [37] neglect Tq Q because they associate it 
with the energy density of the scalar field, which they argue must be uniform 
throughout the Solar System, not larger than the mean cosmological energy 
density and negligible relative to p. Moreover, from Eq. (TT5l) . we see that Tq 
is at least quadratic in $ and can thus be neglected, provided 3/9/(2p/x 2 ) 1. 

Let us now specialize the above treatment to that of homogenous rotating 
source, with mass current 

J = p[ujXr\Q(R-r), (146) 

where R is the radius of the rotating body, p is its density, to is its angu- 
lar velocity, r is the distance from the origin, and is the Heaviside step 
function. One can now use this mass current in Ampere's law to solve this 
equation for A 1 , imposing continuity and finiteness at the origin. The actual 
expressions for the four- vector potential can be found in [37], but its associ- 
ated gravitomagnetic field is given by B = Bqr + -£>cs, where the GR piece is 
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simply 



^ _ AvGptf [ (5 - 3g) h? + 3^f x (fx J), r<i?, 
GR 15 1 ^ [2cJ + 3f x (r x £)] , r > i?, 1 J 



while the CS correction is given by 



B cs = AitGpR 2 {-Di(r) w + D 2 (r) r x + L> 3 (r) fx(fxw)}. (148) 



The function -Di,2,3 (?") are actually functionals of spherical Bessel functions of 
the first je(x) and second yi(x) kind, given explicitly in [37]. Remarkably, the 
gravitoelectromagnetic analogy allows for a solution to the linearized modi- 
fied field equations for an extended rotating source, which possesses both a 
poloidal and toroidal gravitomagnetic field. In other words, the metric does 
not only contain gravitomagnetic components co-aligned with the spin axis 
of the rotating extended body, but also along other axis, perpendicular to 
the plane defined by to and r. Such terms cannot indeed be removed via a 
coordinate transformation since is generically non-vanishing. 

The solution for the metric perturbation found in the gravitoelectro- 
magnetic analogy for extended sources differs slightly from that found in 
the PN scheme for point particles. Some of these differences arise because 
in the latter the point particles were allowed to possess non-vanishing angu- 
lar and linear momentum. The main difference, however, is not due to the 
stress-energy tensor studied, but to the fact that extended and point particle 
treatments are not equivalent in CS modified gravity. In the former, an ad- 
ditional oscillatory behavior is needed, encoded in the appearance of Bessel 
functions, in order to guarantee continuity across the surface of the sphere. 
Such terms are necessary because the C-tensor modifies the junction condi- 

The gravitoelectromagnetic analogy also sheds some light on the study of 
exact solutions to CS modified gravity, reviewed in Sec. HI The CS correction 
to Ampere's law, Eq. (11441) changes the character of the differential system, 
from a first-order one to a second-order one. Such a result is in fact expected 
from the structure of the modified field equations [Eq. (TTBl) ]. since the C- 
tensor depends on third-order derivatives of the metric. One then could also 
expect that approximate solution to CS modified theory need not be easily 
obtained as deformations of GR solutions. This is because, in principle, the CS 
correction could produce dramatic changes to the dynamical behavior of the 
solution (ie. the solution, for example, could be an element of CS \ 7 in the 
nomenclature of Sec. HI which would hold little resemblance to GR solutions). 
CS modified theory and GR, however, have both now been sufficiently tested 
to allow for a perturbative treatment of the modified theory, at least in the 
Solar System. Whether such treatment is allowed in more non-linear, strong 
field scenarios remains to be studied further. 
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5.4 Perturbations of the Schwarzschild Spacetime 

Black hole perturbation theory has been incredibly important in GR, 
leading for example to a better understanding of subtle issues related to BH 
physics and the extraction of gravitational waves in dynamic spacetimes. More- 
over, through the discovery of quasinormal ringing, new tests of GR have been 
proposed (see eg. |124] for a review). These accomplishments alone suggest a 
study of BH perturbation theory in CS modified gravity could lead to inter- 
esting results. Such a study was carried out by Yunes and Sopuerta [M] in the 
non-dynamical formalism with a = k and (3 = 0. 

Let us begin by introducing the basics of BH perturbation theory, as 
developed in GR. Consider then the following perturbed metric 

9ab = 9ab + Kb , (149) 

where g^b is the background Schwarzschild metric and h a b is a generic metric 
perturbation. Henceforth, objects associated with the background will be de- 
noted with an overhead bar. The spherical symmetry of the background allows 
one to expand the metric perturbation in tensor spherical harmonics, thus sep- 
arating the angular i = {9, 0} and temporal-radial /i = {t, r} dependence in 
the perturbed field equations and yielding a 1 + 1 PDE system. This system 
can be further simplified by distinguishing between polar and axial parity har- 
monics: polar (axial) modes acquire a (—1)' [(— 1)' +1 ] prefactor upon parity 
transformations. The aforementioned simplification comes about because in 
GR the perturbative field equations decouple into two subsystem that can be 
classified by their parity. 

With these considerations in hand, we split the metric perturbation into 
h a b = h^ b + h^ b , where each of these pieces is decomposed in tensor spherical 
harmonics: the temporal-radial sector becomes 

h% = 0, = E tft" = E Y £m , (150) 

i,m £,m 

the temporal-angular sector is given by 

h% = E h T = E hi m s e r, h* = e = E?i m >f m , (151) 

£,m i,m i,m £,m 

and the angular-angular sector becomes 

h % = E h T = E Him slT , h l = E h f m = E r2 { K£mY t + G em zfp) , 

i,m £,m £,m £,m 

(152) 

asterisks denote components given by symmetry and where Y lm are standard 
scalar spherical harmonics [see [125] for conventions], Yf™ 1 and S , f m are polar 
and axial vector spherical harmonics respectively, and Y^ m , Zf™ and S^™ 1 are 
polar, polar and axial, tensor spherical harmonics. Vector spherical harmonics 
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are defined for i > 1 via 

yim _ yfen ^ _ y&n _ (153) 

while tensor spherical harmonics are defined for I > 2 via 

yfyn . = y ? ^gn . = y£rn + ^+l) y^.. ? ^gn . = ^ ? ^ 

where Qjj and e 2 - 7 are the metric and Levi-Civita tensor on the 2-sphere re- 
spectively, while colon stands for covariant differentiation on the 2-sphere. All 
metric perturbations are functions of t and r only. 

The perturbative field equations can be decoupled in GR in terms of 
complex master functions, known as the Cunningham-Price-Moncrief (CPM) 
master function [126J and the Zerilli-Moncrief (ZM) master function 
The former is given by 



V CPM — ^ \ n r,t 'H,r 



+ -h e t m ),, (155) 



where = (£ + 2)(£ — l)/2, while the latter is given by 



*zm = W m + (i + h)G im + £ 



fh e ™-rK e r m --(l + \ e )p i r m 
r 



(156) 

where / = 1 — 2M/r is the Schwarzschild factor and Ai = X e + 3M/r. 

The perturbative field equations can be written in terms of these master 
function, obtaining the so-called master equations 



-dl + d% - Vr^\r)\ ^ M/ZM = , (157) 

where is the tortoise coordinate r + = r + 2Mln [r / (2M) — 1] and y f Polar / Axial 
are potentials, which depend on r, parity and harmonic number [125]. In 
principle, the right-hand side of the master equations is not zero, but depends 
on some matter sources, which we neglect here since we are searching for 
vacuum metric perturbations. Moreover, although the master functions are 
gauge invariant, the analysis is simplified if one picks a gauge. Henceforth we 
choose the Regge- Wheeler gauge in which H £m = = G £m = pf m . 

With all this machinery, we can now study BH perturbation theory in CS 
modified gravity. Yunes and Sopuerta |125j studied the non-dynamical for- 
malism, in which the Pontryagin plays a critical role. In terms of perturbation 
functions this constraint reduces to 



*RR 



96M 



vim I / -Lira vim 

'H r, \' l r,t a t,r 



E(£+l)Y £m . (158) 



Although this constraint is automatically satisfied by the background, it is 
not satisfied for generic metric perturbations. Remarkably, this precise com- 
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bination of metric perturbations can be written in terms of the CPM function 
exactly: 

*RR= _ 24M (l+2)^ m y tm (15Q) 

r e (£-2)! CPM ' y ' 

The Pontryagin constraint then forces the CPM function to vanish for all har- 
monics £ > 2, which is not equivalent to requiring that all axial perturbations 
also vanish. One then arrives at the constraint 

h% = h$-*hf». (160) 

The set of allowed solutions is thus reduced by this constraint, which in fact 
tends to lead to an overconstrained system of perturbative field equations, as 
we shall see next. It is also in this sense that the non-dynamical framework 
leads to an overconstrained system of partial differential equations. 

Let us now concentrate on the perturbative field equations in CS modified 
gravity. In order to simplify calculations, Yunes and Sopuerta [125J chose a 
scalar field of the form d = $(t, r), which possesses the same symmetries as the 
background. After harmonically decomposing the perturbative field equations, 
one finds [125j 
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where t/p™ ar and U ^ igX denote polar and axial metric perturbations respec- 
tively 

= , , K £m , G £m ) , C7t ial = (h^ , H em ) , (164) 

and where the right hand sides of Eq. f)16ip and (I163p also depend on deriva- 
tives of The functionals in Eqs. (I16ip and f)163p are long and unilluminating, 
so we shall omit them here, but they are presented in full detail in [125] . Per- 
haps not too surprisingly, the parity- violation induced by CS modified gravity 
breaks the axial-polar decoupling expected in GR. Instead, we now find that 
modes with opposite parity are coupled and cannot in general be treated sep- 
arately. 

Due to the non-decoupling of the perturbative field equations, a generic 
study of their solution is a quixotic task. Yunes and Sopuerta |125j investigated 
several specific cases, from which one can then extrapolate generic conclusions. 
Let us then first consider the canonical choice of d and single-handed metric 
perturbations, ie. purely polar or purely axial perturbations. In either of these 
cases, one can show that the Pontryagin constraint leads to an overconstrained 
PDE system, and thus, if one set of perturbations vanishes, then all metric 
perturbations must vanish. Moreover, these conclusions do not only hold for 
the canonical choice of d but also for many other members of the family 
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$(t,r) [125] . However, a generic result for arbitrary d or arbitrary axial and 
polar metric perturbations has not yet been produced, due to the incredible 
complexity of the perturbative equations. 

The overconstraints of the non-dynamical formalism are mainly due to 
the Pontryagin condition, which is relaxed in the dynamical formalism. Yunes 
and Sopuerta [125] have studied BH perturbation theory in the dynamical for- 
malism, where d and h a b are treated as independent perturbation parameters. 
Consider then the split 

9ab = 9ab + £ h abi ( 165 ) 

= T ($ + e 8$) = t$ + re £ ¥ m Y lm , (166) 

e>i,m 

where here e and r are bookkeeping, independent perturbation parameters. 
The expansion for $ is a bivariate expansion, where $ is some background 
value for the CS coupling and is spherically symmetric, while •& are harmonic 
coefficients of the d perturbation. 

In the dynamical formalism, the Pontryagin constraint is replaced by 
an evolution equation for the scalar field, which to leading order is simply 

= 0, where we have set the potential to zero. To first order in e, the 
evolution equations become 

-ej6{*RR) =er/3 jofr?- 

(167) 

where 5(*RR) is the functional coefficient of *RR to 0(e) given explicitly 
in [125] . Similarly, to leading order in e the perturbative field equations also 
become modified, with corrections that arise from the expansion of the stress- 
energy tensor of By relaxing the Pontryagin constraint, the PDE system 
ceases to be overconstrained and might allow generic metric oscillations. How- 
ever, since rCl, the magnitude of the CPM master function will be con- 
strained to be small, which could lead to interesting observational conse- 
quences, for example in the emission of energy by gravitational wave radi- 
ation. A full dynamical study of the solution to these equations has not yet 
been carried out. 

5.5 Slowly Rotating Kerr-like Black Holes 

Instead of studying arbitrary perturbations, let us now concentrate on 
perturbations that represent slow-rotation deformations of the Schwarzschild 
spacetime. By searching for such solutions, we might gain some insight on how 
to extend the Kerr metric in CS modified gravity. In the previous section, we 
saw that a CS scalar of the form d = $(t,r) will lead to an overconstrained 
system of equations, so in this section we shall explore more general choices. 

The slow-rotation limit of stationary and axisymmetric line elements in 
non-dynamical CS modified gravity (a = k and (3 = 0) was first studied 
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by Konno, et. al. [98]. Consider then the following general, stationary and 
axisymmetric line element 

ds 2 = -f [1 + h(r, 6)} + f- 1 [1 + m(r, 9)\ 

+ r 2 [1 + fc(r, 6)} [dO 2 + sin 2 - w(r, #)dt) 2 ] (168) 

where / = 1 — 2M/r is the Schwarzschild factor, h,m,k,uj are unknown 
functions of r and 0. The metric perturbations are assumed to be linear in 
the perturbation parameter: J/M 2 , the ratio of the angular momentum of the 
rotating compact object to the squared of its mass. With this line element, 
the Pontryagin constraint leads to a condition on the function lu, namely 
u >r g + 2 cot 9uj iT = 0, which leads to the solution to — ui{r)J sin 2 6. Such a result 
is so far independent of the choice of CS scalar. 

The modified field equations can be now linearized in the unknown func- 
tions and solved, given some choice for When a canonical scalar field is 
chosen, Konno, et. al. [98J showed that the linearized equations force these 
functions to vanish and thus a rotating solution cannot be found to first or- 
der. Note that this result is not in disagreement with the discussion in Sec. 15.21 
and 15.31 since there the far- field solutions found cannot be put in the form of 
Eq. (HSU). 

A solution to the linearized CS modified field equation can in fact be 
obtained in the slow-rotation limit, provided we explore other choices for d. 
Konno, et. al. 198 1 showed that with the scalar 



tf = rcos#/A , (169) 

a solution can indeed be found and it is given by 

h(r,9)=m{r,6) = k(r,6) = 0, 
qM = ^/ + ^ \r 2 - 2Mr - AM 2 + AMrf ln(r - 2M)] (170) 

which leads to the Schwarzschild metric, plus a new term in the t<p sector of 
the metric, namely, 

gt , = D 3 f + — \r 2 - 2Mr - AM 2 + AMrf m(r - 2M)] , (171) 

where Di are constants of integration that are assumed linear in J. Note, how- 
ever, that if 7^ 0, the gravitomagnetic sector of the metric naively looks as if 
it could diverge in the limit as r — > 00 increases. On closer inspection, however, 
one finds that invariants and physically relevant observables do not diverge. 
For example, the scalar invariant R abccd R a b cc i oc 48M 2 /r 6 — 4_D 2 /(r 4 sin 2 9), 
which indeed vanishes at spatial infinity, where the divergence at 9 = or it 
presumably arises due to the first-order linear perturbation scheme |129j . The 
quantity Aq in $ is a constant with units of inverse length, which curiously 
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does not appear in the solution for Q. This is because the embedding coordi- 
nate can be factored out and does not enter into the linearized modified field 
equations to leading order in the angular momentum. 

Interestingly, the above solution cannot be interpreted as a small defor- 
mation of the Kerr line element. That is, there is no choice of Di for which g t $ 
can be considered a small deformable correction to Kerr. Such an observation 
implies that the frame-dragging induced by such a metric will be drastically 
different from that predicted by the Kerr line element, in fact sufficiently 
so to allow for an explanation of the anomalous velocity rotation curves of 
galaxies which we shall discuss further in Sec. [71 At the same time, however, 
Solar System experiments have already measured certain precessional effects 
in agreement with the GR prediction [130|I131|I132| . and thus a drastically dif- 
ferent frame-dragging prediction might be in contradiction with these Solar 
System tests. 

Recently, Yunes and Pretorius [51] have extended and generalized this 
result. They showed that in fact the solution in Eq. (11701) is preserved for any 
$ in the family 

$s en = A + A x r cos sin 9 + A y r sin sin 9 + A z r cos 9, (172) 

where are constants. In fact, we can rewrite this CS coupling field as 
$ = babA a x b , where x a = [l,x,y,z] and 5 a b is the Euclidean metric. Note, 
however, that the stress-energy tensor associated with any member of this 
family [including Eq. (11691) ] is constant, and thus the energy associated with 
such a field is infinite. Because of this, the solution found here cannot be 
extended to the dynamical framework. 

Moreover, Yunes and Pretorius [51] also found another solution to the 
slow-rotation limit of the modified field equations, if one considers a generic 
CS scalar field in the non-dynamical framework: 

= f( r , 0) + rg((j)) + rh (d<f> - t) + rk(9, 0) 
dr 



_l_ r 



r 



-d r f(r, 0) + -f(r, 0) + ^j{r) 



where /, g, h, j and k are arbitrary functions and C\ is another integration 
constant. This arbitrary function can be chosen such that the new CS scalar 
possesses a sufficiently fast decaying stress-energy tensor with non-infinite 
energy. For example, iff — g — h — k — and j = —3jo/r 2 , then •& = 
jo/r 2 , for constant j , and thus Eq. (11731) is compatible with the dynamical 
framework. 

The existence of two independent solutions to the modified field equations 
in the non-dynamical framework and in the slow-rotation limit suggests that 
there is a certain non-uniqueness in the framework encoded in the arbitrari- 
ness of the choice of This scenario is to be contrasted with the dynamical 
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framework, where the CS scalar is uniquely determined by its evolution equa- 
tion and there is no additional freedom (except for that encoded in initial 
conditions). 

In view of this problems with the non-dynamical framework, Yunes and 
Pretorius [M] studied the same scenario but in the full dynamical framework. 
A new approximation scheme is employed on top of the slow-rotation require- 
ment, which essentially demands that the CS correction be a small deformation 
of the Kerr line element, ie. the CS coupling, is assumed small relative to the 
GR one. In this way, one finds the solution_J [84J: 



where ds^ eTT is the slow-rotation limit of the Kerr metric, M is the BH mass 
and J = a ■ M is the BH angular momentum to leading order. This solution 
is valid to second order in the slowly-rotation expansion parameter a/M, as 
well as in the strength of the coupling a/ (I3kM a ). Notice that this solution, 
derived under the dynamical formulation, is perfectly well-behaved at spatial 
infinity, remaining asymptotically flat. 

Equation (11741) is the first0 rotating BH solution in dynamical CS modi- 
fied gravity, and can be thought of as a small deformation of a Kerr black hole 
with additional CS scalar "hair" of finite energy. Although this is a "hairy" 
solution, Sopuerta and Yunes [99| have shown that the solution is still entirely 
described by the mass and the angular momentum of the source. The no-hair 
theorem is, however, violated in that the relation between higher-multipoles 
and the mass quadrupole and current dipole is CS modified at t = 4 multipole 
due to the CS correction in Eq. (I174p of the gravitomagnetic sector. 

The solution found in the dynamical theory presents interesting parity 
properties. Since the dynamics of the CS scalar field are determined by the 
Pontryagin density, this field is parity- violating (ie. it is a pseudo-scalar). Both 
the Kerr metric and the CS correction to it are also parity-violating, but the 
later is induced by a curvature-scalar field interaction, instead of due to the 
Kerr distributional stress-energy. In fact, the parity violation introduced by 
the CS correction becomes dominant in regions of high curvature. 

Moreover, the dynamical solution presented above shows remarkable sim- 
ilarities with some of the far-field results found in the non-dynamical frame- 



7 Notice that the solution for $ is identical to that found by Campbell |85| and 
Reuter [57] and discussed in Eq. ([90]) . except that there the backreaction of this 
field on the metric was ignored. 

8 Shortly after publication of this result, Konno, et. al. employed slightly different 
methods to verify that the solution found by Yunes and Pretorius indeed satisfies 
the modified field equations |133j . 




(174) 
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work. In particular, it is only the gravitomagnetic sector of the metric that is 
CS modified, thus implying that the frame-dragging effect will be primarily 
corrected. Notice, however, that the CS correction is highly suppressed by large 
inverse powers of radius, which suggest that weak-field tests will not be able 
to constrain the dynamical framework. We shall investigate this possibilities 
further in Sec. [71 



5.6 Gravitational Wave Propagation 

GW solutions have been studied by a large number of authors [221251251331351511591 , 
but mostly in the non-dynamical theory, which we shall concentrate on here. 
The first GW investigation in non-dynamical CS modified gravity was car- 
ried out by Lue, Wang and Kamionkowski [134] . who studied the effect of 
GWs in the cosmic microwave background (see Sec. [5J). Jackiw and Pi [22] 
also studied GWs in CS modified gravity, concentrating on the generation of 
such waves and the power carried by them in the modified theory. Shortly 
after, such waves were used to explain baryogenesis during inflation [28J and 
to calculate the super-Hubble power spectrum [29]. The generation of GWs 
was also studied in the dynamical formalism [81] through the construction 
of an effective stress-energy tensor and the Isaacson scheme [95] • Recently, 
GW tests have been proposed to constrain CS gravity with space-borne [33J 
gravitational wave interferometers. 

We shall here discuss GW solutions in non-dynamical CS modified gravity 
and postpone any discussion of GW generation to the next section. Moreover, 
we shall not discuss in this section cosmological power spectra, since these 
will be summarized in Sec. |SJ GW propagation in CS gravity has only been 
studied in the non-dynamical formalism with (3 — and a = k. Let us begin 
with a discussion of GW propagation in a FRW background. Consider then 
the background 

ds 2 = a 2 (r]) [-air] 2 + + d^d^] , (175) 

where a(rj) is the conformal factor, rj is conformal time and \ l are comoving 
coordinates. The quantity stands for the gravitational wave perturbation, 
which we take to be transverse and traceless (TT), h := h l i = S^hij = and 
dih 1 ! = 0. One can show that a coordinate system exists, such that the gravi- 
tational wave perturbation can be put in such a TT form. For the remainder 
of this section, i, j, k stand for spatial indices only. 

With such a metric decomposition, one can linearize the action to find 
the perturbed field equations. In doing so, one must choose a functional form 
for the CS scalar and we shall here follow Alexander and Martin [21], who 
chose $ = $(•)]). One can show that the linearized action (the Einstein-Hilbert 
piece plus the CS piece) to second order in the metric perturbation yields 
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Seh + S C s = t / d A % 
4 jv 



k 



a 2 (r]) (h l j^ h J itV - h l j>k 
- 0,„ l%3k (h q i>v h kq , jv - h%, r h kq , rj )] + 0(h) 3 (176) 



where a — k, commas in index lists stand for partial differentiation, rj is 
conformal time and e 1 ^ = e TOfc . Variation of the linearized action with respect 
to the metric perturbation yields the linearized field equations, namely [29J 

UhU := ie**G [(<d m - 2m tV ) h l)Km + # v nh i)kjP ] , (177) 

where □ is the D'Alembertian operator associated with the background, namely 

i :/ /. w • •_>:)(/., - (178) 

with / some function of all coordinates and the conformal Hubble parameter 
"K := a )V /a. One could have, of course, obtained the same linearized field 
equations by perturbatively expanding the C-tensor. 

One can see from Eq. (I177P that the evolution of GW perturbations is 
governed by second and third derivatives of the GW tensor. Jackiw and Pi [22] 
were the first to point out that for the canonical choice of $ the GW evolution is 
governed by the D'Alembertian of flat space only, if we neglect corrections due 
to the expansion history of the Universe (ie. if this vanishes, then the linearized 
modified field equations for the GW perturbation are satisfied to linear order). 
Such a result implies there are two linearly independent polarizations that 
propagate at the speed of light. 

Let us now concentrate on gravitational wave perturbations, for which 
one can make the ansatz 

A - ■ 

hij = -^r exp [-i [cj)(r]) - Kn kX k )\ , (179) 

where the amplitude Aij, the unit vector in the direction of wave propagation 
n k and the conformal wavenumber n > are all constant. It is convenient to 
decompose the amplitude into definite parity states, such as 

= A R e% + A L e% (180) 

where the circular polarization tensors ef- L are given in terms of the linear 
ones ep x by [21] 
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These polarization tensors satisfy the condition 



n i e«*ej5 L = iA R ^(e'',) R,L , (183) 



where Xr = +1 and Al = — 1 

With this decomposition, the linearized modified field equations [Eq. (11771) 
reduce to 



a 2 



(184) 



Before attempting to solve this equation for arbitrary $(77) and a(rj), it is 
instructive to take the flat-space limit, that is a — > 1, and thus, a — ► 0. 
Assuming further that time derivatives of the CS scalar do not scale with !K, 
one finds that the above equation reduces to [2"91I3"3"] 

(^r,l + 0r, l - k 2 ) (l - A R , L A;tf) = zA R)L A; RL , (185) 

where k is the physical wavenumber 3-vector, t stands for cosmic time and 
overhead dots stand for partial differentiation with respect to time. Let us 
further assume that the GW phase satisfies jW /0 2 t <C 1, which is the standard 
short- wavelength approximation, as well as •& = $q = const. Then the above 
equation can be solved to first order in •& to find 

cfyit) = O + kt + l -^^kt + O(tf) 2 , (186) 

where <po is a constant phase offset and the uncontrolled remainder 0(i?) 2 
stands for terms of the form k 2 d 2 or kdd. The imaginary correction to the 
phase then implies an exponential enhancement /suppression effect of the GW 
amplitude, as this propagates in CS modified gravity. Recall that here we are 
interested in the propagation of GWs, which is why the right-hand side of 
Eq. (1185P implicitly omits the stress-energy tensor. We shall see in the next 
section that if there is a stress-energy tensor, then the CS correction depends 
both on the first and second derivatives of the CS scalar. Lastly, if we had not 
assumed that d = const, then the solution would have become 

0(t) = ( p + kt + l -^±M(t) + 0(tf) 2 , (187) 

which still exhibits the exponential suppression/enhancement effect. 

The exponentially growing modes could be associated with instabilities in 
the solution to gravitational wave propagation. One must be careful, however, 
to realize that the results above have been obtained within the approximation 
k 2 $ 2 < 1 > k$$. Thus, provided i? is smaller than the age of the universe, 
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then the instability time scale will not have enough time to set in. For larger 
values of the approximate solutions we presented above break down and 
one must account for higher order corrections. A final caveat to keep in mind 
is that these results are derived within the non-dynamical formulation of the 
theory; GW solutions in the fully dynamical theory are only now being actively 
investigated. 

Let us now return to the field equation for the phase in an FRW back- 
ground [Eq. f)184p ]. The solution to this equation is now complicated by 
the fact that the scale factor also depends on conformal time, and thus, 
one cannot find a closed form solution prior to specifying the evolution of 
a(rj). Let us then choose a matter-dominated cosmological model, in which 
a (77) = a Q T] 2 = ao/(l + z), where ao is the value of the scale factor today and z 
is the redshift. It then follows that the conformal Hubble parameter is simply 
given by % = 2/r] = 2(1 + z) 1 / 2 . With this choice, one can now compute the 
CS correction to the accumulated phase as the plane-wave propagates from 
some initial conformal time to rj, namely 



A0(r,l) = i\R,LkH [ 



(188) 



which one can check reduces to the flat space result of above in the right limit. 
Note that the exponential enhancement /suppression effect now depends on an 
integrated measure of the evolution of the CS scalar and the scale factor. 

The CS correction to the GW amplitude derives from a modification to 
the evolution equations of the gravitational perturbation, but it also leads 
to important observational consequences. One of these can be understood by 
considering a GW generated by a binary black hole system in the early inspiral 
phase. The GW produced by such a system can be described as follows 



M 



2/3 



1 + Ar, l cos if exp [-i (tf (t) + A0 RjL ) 



;i89) 



where d^ = a^r] (1 + z) is the luminosity distance to the binary's center of 
mass, is the GW phase described by GR, k (t) is the instantaneous wave 
number of the gravitational wavefront passing the detector and 3YC is the co- 
moving chirp mass, which is a certain combination of the binary mass com- 
ponents. The inclination angle 1, the angle subtended by the orbital angular 
momentum and the observer's line of sight, can be isolated as 



1 + COS I 



COS I 



exp 



2k(t) 



c 



1 + COS I 

1 — cost' 



(190) 



where we have defined 



Hi 



drj 

7?' 



(191) 
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We see then that from an GW observational standpoint, the CS correction 
leads to an apparent inclination angle I, which effectively modifies the actual 
inclination angle by a factor that depends on the integrated history of the CS 
correction: 

_ sinh(M|k)) +cosh (MgM) COS6 UM(zl . 2 , , 

cos t = >, .'. / N jt ', .( ~ cos t-\ — sir l+0 4 . 

cosh (M|M) + sinh (^fM) cos i Ho K J 

(192) 

We see then that the CS correction effectively introduces an apparent evolution 
of the inclination angle, which tracks the gravitational wave frequency. 

The interpretation of the CS correction as inducing an effective inclina- 
tion angle should be interpreted with care. In GR, if a gravitational wave 
propagates along the line of sight, such that the actual inclination angle is 
(eg. or 7r), then the amplitude is a maximum. In CS gravity, however, the 
amplitude can be either enhanced or suppressed, depending on whether the 
wave is right- or left-circularly polarized. When the CS effect suppresses the 
GW amplitude, one can think of this as an effective modification of the incli- 
nation angle away from the maximum. However, when the CS effect enhances 
the amplitude, there is no real inclination angle that can mimic this effect 
(ie. the effective angle would have to be imaginary). 

The evolution equation for the gravitational wave perturbation depends 
sensitively on the scale factor evolution [see, eg. Eq. (j!77j) ] . Alexander and 
Martin [29] have investigated gravitational wave solutions when the scale fac- 
tor presents an inflationary behavior. Suffice it to say in this section that 
Eq. ( 1177ft can be recast in the form of a parametric oscillator equation, with 
a non-trivial effective potential. In certain limits appropriate to inflation, one 
can solve this differential equation in terms of Whittaker functions, which can 
be decomposed into products of trigonometric functions and exponentials. In 
essence, the solutions present the same structure as that of a matter-dominated 
cosmology. Once the gravitational wave modes have been computed, one can 
proceed to calculate the power spectrum, but these results will be discussed 
further in Sec. 18.11 



5.7 Gravitational Wave Generation 

The issue of GW generation by dynamical matter sources in CS modified 
gravity was first studied by Jackiw and Pi [22] ■ Once more, this problem has 
been treated only in the non-dynamical formalism ((3 = and a = k) and 
with the canonical choice of although very recently the much more difficult 
problem of GW generation in dynamical CS modified gravity has begun to be 
investigated [99] . 

In the presence of a stress-energy tensor, the modified field equations 
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linearized about a Minkowski background (g a b = Tjab + h a b) become 



U v hh + ^U v h i)kyP - $<? k{ % Kp = ~-T v (193) 

Kj 

where t is the standard time coordinate of Minkowski spacetime, e pjfc = e° rafc , 
the D'Alembertian operator is □ := — df + d^didj and overhead dots stand 
for partial differentiation with respect to time. One can derive this equation 
simply from Eq. (I177p by taking the limit a — > 1 and !K — > and reinserting 
the stress-energy tensor T^, which must now be TT, since so is hij implicitly. 

Although the explicit solution to the above equation has not yet been 
computed for a general CS scalar field, this problem has been studied for the 
canonical choice of With the assumption that = 0, we then obtain, to 
first order, the solution found in Eq. (11201) . which for a TT stress energy can 
be recast as 

Dnhii = (194) 

K 

where the effective stress-energy has been defined as 

% := Tij - re k \{T j%k . (195) 

We see then that GWs generation in non-dynamical CS gravity with a canoni- 
cal d is nothing but GR GWs in the presence of such an effective stress-energy 
tensor [22] • 

For concreteness, let us assume the stress-energy tensor represents a 
monochromatic source, with definite frequency u, which is radiating GWs 
in the z-axis with wave-vector k. The only non-vanishing components of the 
stress-energy tensor are then T xx = —T yy and T xy = T yx and it's Fourier trans- 
form shall be denoted T + and T x respectively. We then find that the effective 
stress-energy becomes 

fij (T + -ikiT x T X+ «T + 
\T x +iMT+ -T++iMT x 

which exhibits the natural mixture of polarizations of CS GWs. 

A natural next step is to compute the power carried by such CS modified 
GWs per unit angle dVt. Since CS GW theory for a canonical d is identically 
equivalent to GR with an effective stress-energy tensor, it follows that the GW 
power is given by 

HP 

= 16kG 2 uj 2 f *ijf ij . (197) 

One is tempted to insert the solution for the effective stress-energy found in 
Eq. (I196P into Eq. (11971) . This effective stress-energy, however, is insufficient, 
since it is only a solution to first-order in t? and the expression for the power 
emitted requires a second-order solution. Jackiw and Pi [22J have found this 
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solution, which reduces to Eq. (11 96ft multiplied by {l — k 2 d 2 ) 1 . One then finds 
that the power emitted is given by 



dR 



R,L 



dn 



32kG V 



k 2 $ 2 



+ 2iM (t + T* -T X T] 



1 + k 2 d 2 ) (T 2 + T. 



(198) 



which when linearized to second order in d becomes 



dP» 



^ ~ 32kG 2 u 2 [(l - k 2 $ 2 ) (T 2 + T 2 ) + 2z&tf (T+T* - T X T* + )] . (199) 



The power carried by circularly polarized GWs (T+ = iT x ) is corrected by CS 
gravity only to second order in d [22J. 

The results presented above hold only for the canonical choice of d. Had 
we allowed the second derivative of the CS scalar to be non- vanishing, we 
would have found a linear correction to the power carried by CS GWs. In 
the dynamical formalism, this effect is more clear as the CS scalar also carries 
energy-momentun away from the system [99] . Such modifications to radiation- 
reaction would affect the inspiral and merger of binaries, as suggested in |99j . 
which leads to a powerful test of the dynamical formulation that we shall 
discuss in Sec. [3 



6 Non-Vaccum Solutions and Fermionic Interactions 

The first-order formulation of CS modified gravity was first studied by 
Cantcheff |135j . who realized that the modified theory would lead to non- 
vanishing torsion for a canonical CS scalar. Such an idea was later generalized 
to arbitrary 9 in the non-dynamical formalism and the torsion tensor was 
specialized to Earth's gravitational field [6]. Such a formulation naturally al- 
lows for the coupling of fermions to the modified theory, thus permitting the 
study of non- vacuum spacetimes [6]. We shall summarize these results here, 
beginning with a description of the first-order formalism and Einstein-Cartan 
theory, following mainly [r§fl3b'|ll37| . We then continue with a discussion of 
the first-order formulation of CS modified gravity. 

6.1 First- Order Formalism 

Consider a 4-dimensional manifold M with an associated 4-dimensional 
metric g^. At each point on this manifold, let there be a tetrad e^, so that 
the metric can be recast as g a b = e^e^rjjj, where rjjj is the Minkowski metric. 
Internal indices range I, J — (0,1,2,3), just as spacetime indices do, and 
often, we will suppress spacetime ones and the tetrad shall be written e 1 . We 
can raise or lower spacetime and internal indices with the metrics g a b and r\u 
respectively. 
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The introduction and differentiation of internal versus spacetime indices 
is crucial to the first order formalism. Riemannian fields, like the metric tensor, 
exist on the base manifold M and have a finite dimension, but gauge fields 
can be infinite dimensional and so they must exist on a different vector space. 
For example, the tetrad e 1 and the spin connection uj kl are 1-forms on the 
base manifold, while the curvature tensor associated with it, F KL , is a 2-form 
on the base manifold. On the other hand, the fermion field ip is a 0-form on 
the internal vector space. A fiber bundle is defined as the union of the base 
manifold and the internal vector space, with each fiber a different copy of the 
internal vector space. One can think of the Lie group associated with the fiber 
bundle as glueing all fibers and the base manifold together [19]. 

The recovery of spacetime indices is sometimes achieved via the wedge 
product operator, defined via 

(A^B) ab := { P^A [ai ... ap B h ... bq] (200) 

with A and B p- and g-forms respectively. Note that we shall here follow 
the convention that spacetime indices always appear after internal ones. The 
wedge product operator satisfies the following chain rule 

D {u) (AAB) — (D (u) A) AB + {-If A A (D H B) , (201) 

and the following commutativity relation 

A A B = (—l) pq B A A. (202) 

Also note that since the wedge product acts on spacetime indices only, it exists 
on the base manifold and not on the internal space. 

Now that the metric and tetrad have been defined, let us introduce the 
generalized covariant derivative operator D and the spin connection uj ij . Given 
a tensor A KL a we can define the covariant derivative as 

D H A KL := dA KL + u KM A A M L + u LM A A K M , 

(203) 

D(u>)Akl '■= dA KL — lo k m A A M l — u L M A A KM , 

(204) 

where the exterior derivative operator d acts on spacetime indices only: 

dA KL := 2d [a A KL b] . (205) 

The commutator of covariant derivatives allow us to define the curvature ten- 
sor associate with uj ij , namely 

F KL = duo KL + uj k m A lu ml , (206) 
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which reduces to the Riemann tensor if the spin connection is metric compat- 
ible and torsion-free, ie. if the connection is the Christoffel one. One can then 
show after some calculation that 

5 U1 F IJ = D {u)) 5uo IJ . (207) 

The spin connection has a certain degree of freedom that can be fixed by 
demanding that it be internally metric compatible D^rj IJ = 0. This condition 
forces the connection to be completely antisymmetric on its internal indices 
ijj\ IJ ) = o. We can then define the torsion tensor as 

T 1 := D^e 1 = de 1 + u\t Ae M , (208) 

which is equivalent to T 1 ab = 2D[ a eL. If we reinstate spacetime indices, we 
find 

T a bc = 2C a [bc] , (209) 
where C a bc is the antisymmetric part of the connection, or contorsion ten- 
sor. Note that internal metric compatibility is not equivalent to a torsion-free 
condition. 

The contorsion tensor can be defined purely in terms of wedge products. 
Let us then decompose the spin connection into a symmetric and tetrad com- 
patible piece T 1 j and an antisymmetric piece C 1 j, the contorsion tensor. The 
torsion tensor is then given in terms of the contorsion via 

T 1 = C J j A e J , (210) 

which can be inverted to find 

Cuk = ~\ (T IJK + T JKI + T KJI ) . (211) 

The contorsion tensor is fully antisymmetric on its first two indices, while the 
torsion tensor is fully antisymmetric on its last two indices. Equation (12091) 
can also be obtained from Eq. (12081) in spacetime indices, if we use the trans- 
formation law from spin to spacetime connection: 

u) ab e = ej ujK a e b - ej d a e{, (212) 

which can be established from D^e 1 = 0. Eq. f 1 2 1 2 j) is sometimes referred to 
as "the tetrad postulate" . 

The curvature tensor can be expressed purely in terms of the Riemann 
tensor (that depends only on T 1 j) and terms that depend on the contorsion 
tensor: 

F IJ = R IJ + D {r) C IJ + C I M AC MJ , (213) 
where D^) is the connection compatible with the symmetric connection. The 
Bianchi identities in first-order form become 

D^T 1 = R I K Ae K , D {UJ) R IJ = 0. (214) 
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6.2 First-Order Formulation of CS Modified Gravity 



Let us now apply the formalism of the previous section to GR and to the 
modified theory. The Einstein-Hilbert action can be recast in terms of forms 
as 

Seh = ^ euKLe 1 A e J A F KL . (215) 
2 Jv 

One can convert this into Eq. ([2]) by rewriting the curvature tensor as F IJ = 
(1/2)F IJ KL e K A e L and using the identity 

e 1 A e J A e K A e L = -a e IJKL , (216) 

where a = y/^g d 4 well as the Kronecker-Delta relations 



abed 

£ Cabef O e °/) 



4 6%Sf, (217) 
+A^5^8f, (218) 
V abc % b ef = +A5^5f. (219) 



V^tabef = +4 V-5 Wf, 



Note that full internal index contractions are identically equal to spacetime 
index contractions, if the quantities contracted are tensors. 

Similarly, we can attempt to recast the CS action in terms of forms via 

Scs = 77 / $ F A F, (220) 



v 

where the integrand reduces to RAR = RijAR ij for a symmetric connection. 
Again, Eq. (I220p can be converted into Eq. ([3]) in the same way as above, which 
then establishes that a*RR = 2RAR for a symmetric connection. We can now 
integrate by parts to obtain the CS action in terms of the Pontryagin current, 
but to do so we first must realize that F A F = dQ^, where the Chern-Simons 
3-form is defined as [28 J27|31] 

2 

tt 3 = u IJ A duj u + -co 1 j A uj j K A tu K /. (221) 
3 

We then find that Eq. (12201) can be rewritten as 

Scs = ~ [ *?Afi 3 , (222) 
2 Jv 

where we have neglected the boundary contribution. Converting this expres- 
sion to spacetime indices, one can check that one recovers Eq. ( JTUl) if the 
connection is symmetric. In fact, the dual to this 3-form is actually the Pon- 
tryagin current in Eq. ([9]) in disguise, while *(dQs) = *FF/2. 

We can now vary the action with respect to all degrees of freedom to 
obtain the modified field equations. The variation of the action with respect 
to the tetrad yields 

1 2 

Fab — 7^9abF = —T ab , (223) 

Z K 
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since the CS current does not depend on the tetrad, but only the spin con- 
nection. Equation (I223P resembles the Einstein equations, except that the 
quantity F ab = F c acb and F = F a a are not the Ricci tensor and scalar, but 
contractions of the curvature tensor. 

The full curvature tensor can be reconstructed once one solves for the 
contorsion tensor via the torsion condition of the modified theory. This con- 
dition is arrived at by requiring that the variation of the action with respect 
to the spin connection vanishes. Let us then rewrite the CS 3-form via 

Q 3 = J j A F J j - -J j A uo J K A u K j. (224) 

such that its variation with respect to ujkl reduces to F KL . The variation of 
the action with respect to u KL reduces to 

6S=-~ f D (euKLe 1 A e J ) 5u KL - ^ [ d# A F KL 5u KL , (225) 

where we have set (3 = for simplicity and we have integrated the first term 
by parts. Cantcheff [135j has shown that SS == is equivalent to the torsion 
constraint 

D (euKLe 1 A e J ) = --d-d A F KL , (226) 

which is nothing but the CS modified second Einstein-Cartan structure equa- 
tion. This equation can also be rewritten in terms of the torsion tensor as 

e IJKL T 1 A e J = -^d$ A F KL . (227) 

The torsion condition of Eq. (I227p shows that in CS modified gravity the 
connection must be torsionfull. Solving for the torsion in Eq. (I227P is non- 
trivial, due to the implicit appearance of torsion itself through the curvature 
tensor on the right-hand side of this equation. However, if we replace this 
tensor by the Riemann (torsionless) tensor, then we can solve for the torsion 
tensor exactly to fincf^l. 

T i = « € ijkl vj Rr ^ (22g) 

where here we have used the first Bianchi identity R\ijk]l = 0. Equation (12281) 
should be interpreted as an approximate solution to the torsion constraint, 
where we have neglected higher than second powers of 9 |135|l6] . 

One can now relate this torsion tensor to the gravitational field of a 
compact binary in the PPN formalism [6J. One can show that this tensor is 
proportional to the product of v a with V x V and V£7, where V 1 and U are 



9 This solution is the same as that found by [135] with a = 2 and shortly after 
by [U] with a = k. We have further checked that Eq. (|228|) is indeed a solution to 
the torsion constraint of Eq. (|227p 
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the PPN vector and Newtonian potentials. In fact, there is no choice of 9 for 
which all components of the torsion tensor vanish. Since this tensor affects 
the motion of point particles (in particular the frame-dragging effect [B]), 
one concludes that non-dynamical CS modified gravity generically leads to 
modified precession, irrespective of 9. Care must be taken, however, since in 
spite of generically being non-zero, the torsion tensor is after all proportional 
to the dual to the Riemann tensor. Thus, Earth-based experiments that search 
for non- vanishing torsion [13 8113911140] cannot measure this effect, since the 
Riemann tensor on Earth is prohibitively small. 

Are the first and second-order formulations of CS modified gravity equiv- 
alent? The modified field equations of Eq. (I223j) resemble the Einstein equa- 
tions, except that the curvature tensor contains torsional pieces. When Eq. (12281) 
is used as the torsion, one can show that all these torsional pieces conspire 
to produce the C-tensor of Eq. (I18p |135j . One then finds that the first and 
second-order formulations of CS gravity are equivalent if and only if the CS 
action in first-order form is defined in terms of the torsion-free curvature ten- 
sor (ie, in terms of the symmetric connection), such that Eq. (12281) is the exact 
solution to the torsional condition. 

The CS 3-form, however, has been here introduced in terms of the gener- 
alized spin connection, for which the torsion condition is non-linear, depending 
explicitly on the curvature tensor. For such a torsion condition, the torsion 
tensor in Eq (I228j) is still formally valid but only to linear order in 9, and thus, 
the first-order formalism is still equivalent to the second-order one but only 
to 0($). The 0(-#) 2 corrections to Eq. (I228j) will modify the field equations, 
and thus, force them to not be equivalent to those of the second-order formal- 
ism. Cantcheff [135J has further shown that line elements that are solutions 
in the second-order formulation of CS gravity (eg. the Schwarzschild metric, 
which leads to a vanishing C-tensor) are not necessarily solutions to the field 
equations of the first-order formulation if higher-order in 9 terms are included. 

The inequivalence between the first- and second-order formalism thus de- 
pends on quadratic or higher powers of the CS scalar. The CS modified action 
considered in Sec. 12.11 however, only considers linear terms in 9. In principle, 
there will be 9 2 and higher-order terms in the action that one could have to 
include, since these will also generically break parity [eg. d 2 ( *RR) 2 }. A con- 
sistent comparison between first- and second-order formalisms thus requires 
that such terms be taken into account, if one wishes to define the CS action 
in first order form in terms of the torsionfull curvature tensor, instead of the 
Riemann tensor. 



6.3 Fermions and CS Modified Gravity 

One of the advantages of the first-order formalism is that it allows for the 
inclusion of fermions and bosons in the action. Let us then considering the 
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inclusion of the following piece to the full action of Sec. 12.11 



Sd = \! d " x ^ (We?^ + c.c.) , (229) 

where c.c. stands for complex conjugation, if) is a Dirac spinor, 7 7 are gamma 
matrices and e is a coupling constant. Fermions are here represented by Dirac 
spinors, which are gauge field that live naturally in SU(2). The tetrad and 
the spin connection of GR, on the other hand, are fields that live in SO(l, 3). 
Thus, without the machinery of fiber bundle theory one could not easily couple 
fermions to CS modified gravity. 

From the group structure of if), one can deduce how the generalized SU (2) 
covariant derivative acts on Dirac spinors: T> a ip := d a ip — (l/4)u IJ alilj if>, 
where we shall here follow the sign conventions of [75]. Variation of the full 
action with respect to Dirac fermions then leads to the massless Dirac equa- 
tion [6]: 

7 °£>i r V = \^ M Ca KL l M lKlL^, (230) 

where notice that we have not included a mass term for the fermions for 
simplicity. We see then that the Dirac equation is modified in CS modified 
gravity by a source term that depends on the contorsion tensor. Equation (12301) 
implies that the CS effects might be enhanced in spacetime regions where the 
momentum of Dirac fermions is large. 

Before we can vary the full action with respect to the connection it is 
convenient to recast it in first-order form. Doing so, Eq. (I229p becomes 

S D = ^j tiJKLe 1 A e J A e K A (# 7 L !D</> + c.c.) . (231) 
Upon variation of the action with respect to the connection, Eq. (12271) becomes 
cijklT 1 A e J = ~dti A F KL + ^-e K A e L A J (5) , (232) 

where J§ := if>'~f5'J L ip is the fermion axial current, e is the determinant of the 
tetrad field and we have used the identity 

7 W 1 = -ie IJKL ^i L + 2^^. (233) 

This equation can be solved for the torsion tensor if we once again replace the 
curvature tensor by the Riemann tensor on the right-hand side. Doing so, we 
find [6] 

T 1 = -^e IJKL vjR KL - fjjKLJfe' A e K , (234) 
4k qk 

which essentially is an approximate solution that neglects quadratic and higher 

powers of 6. From this, the contorsion tensor can be calculated using Eq. ( 12111) 

to find 

Cu = —VNei! NML Rjkjml + -^-e IJKL J^e K . (235) 
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Since the first- and second-order formalism are equivalent to leading order 
in 6, one can compute the interacting action be reinserting the torsion solution 
in the full action. Doing so, one finds [B] 

s = Is + w K L s ( 2J >' fl - - J >° R ) (236) 

where we have neglected other terms that are either higher-order in 9 or in 
the gravitational coupling constant. The first term is the standard 4-fermion 
interaction that arises in Einstein-Cartan theory coupled with fermions. The 
second term is a new CS contribution that depends on the embedding co- 
ordinate, as well as the Ricci tensor and scalar. This new term represents a 
2-fermion interaction and it is not suppressed by the gravitational coupling 
constant. We can then conclude that, at least to linear order, fermion current 
tend to enhance the CS correction, some of the implications of which shall be 
discussed in Sec. [7J 

7 Astrophysical Tests 

All tests of CS modified gravity to date have been performed with as- 
trophysical observations and concern the non-dynamical framework. After 
Alexander and Yunes [7|32] realized that the modified theory predicts an 
anomalous precession effect, Smith, et. al. [37J tested the non-dynamical model 
with canonical CS scalar with LAGEOS [130lll31|ll32j and Gravity Probe 
B |141|142j observations, placing the first, albeit weak, bound on the CS 
scalar. In view of these results, Konno, et. al. [129J proposed that the CS 
correction could be used to explain the flat, rotation curves of galaxies , which 
in turn could yield yet another constraint on the non-dynamical theory for 
non-canonical Recently, Yunes and Spergel [143] used double binary pulsar 
data to place a bound on the non-dynamical model with canonical CS scalar 
that is eleven orders of magnitude stronger than the Solar System one. 

The dynamical model remains untested today, mainly due to the difficulty 
in calculating observable quantities in a consistent way. The only possible 
avenue to perform such a test seems currently to be through gravitational wave 
observations [33]. Cosmological tests of the modified theory will be discussed 
in the next section. 

7. 1 Solar System Tests 

The non-dynamical modified theory has been so far only tested through 
Solar System, frame- dragging experiments. Smith, et. al. [37] studied the 
anomalous precession inherent to the non-dynamical model with the conven- 
tions a = —£/3 and (3 = 1 and we shall summarize these results here. With 
these conventions, d has units of inverse length or mass and mcs is a charac- 
teristic mass scale defined in Eq. ( 11451) . When testing CS modified gravity, we 
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shall in fact place bounds on this parameter, which can be trivially related to 
i? by the inversion of Eq. (11451) . 

Only Solar System tests that sample the gravitomagnetic sector of the 
gravitational field can be used to test CS modified gravity in the non-dynamical 
model. As shown in Sec. |5l the non-dynamical modified theory possesses the 
same PPN parameters as GR, except for the gravitomagnetic potential. This 
implies, in particular, that the perihelion shift of Mercury or light deflection 
by the Sun cannot be used to constrain the modified theory. The only physi- 
cal effect of the non-dynamical CS modification is the induction of anomalous 
precession effects. 

Precession is a generic term used to address the change in the rotation 
3-vector of some spinning object, ie. a non-vanishing 4-gradient of the spin 
angular momentum. Two types of precession can be distinguished: torque-free 
and torque-induced. The former corresponds to situations in which the spin 
angular momentum is not coaligned with the axial Killing vector. The latter, 
also known as gyroscopic precession, occurs in situations where there is an 
additional torque (such as that of a gyroscope) that pushes on the spin an- 
gular momentum vector, forcing it to wobble. Gyroscopic precession can be 
studied in a Newtonian framework, but relativity adds three additional correc- 
tions: Thomas precession, an additional special relativistic correction due to 
the observer's non-inertial rotating frame; de Sitter or geodetic precession, a 
GR effect that accounts for Schwarzschild-like deviations from flat spacetime; 
Lense-Thirring precession, a GR correction due to the gravitomagnetic sector 
of the Kerr metric. 

The CS modification can correct several different types of precession, de- 
pending on the physical scenario under consideration. Torque-free precession 
can occur if one considers the far-field expansion of a non-dynamical CS spin- 
ning black hole metric, where the axis of rotation seems not to be co-aligned 
with the axial axis of symmetry. In this one-body scenario, in the absence 
of external torques, the spin angular momentum of the black hole precesses 
around the symmetry axis in a wobbling fashion. The evolution of the wob- 
ble angle requires the determination of the Killing axis, from which one can 
obtain the frequency of precession, via Uf ?a J/I\ to Newtonian order, where 
I\ is the moment of inertia about the symmetry axis [144] . Since such an ar- 
rangement is asymmetric, one expects GW emission leading to spin-down and 
alignment. In the non-dynamical formalism, then, spinning black holes would 
tend to "unspin" themselves via interactions with the CS scalar, thus relax- 
ing to the Schwarzschild solution. If so, observations of spinning black holes 
could be used to constrain the magnitude of the canonical scalar [145] . Such 
a possibility has not yet been studied in detail. 

Torque-induced precession is also modified by the CS correction through 
the correction to the gravitomagnetic sector of the metric. Consider first the 
motion of a test body in the external field of a CS spinning source. The far-field 
solution for such a source was summarized in Sec.[5]for extended bodies [37]. In 
such a field, the orbital elements of a test body will experience Lense-Thirring 
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precession |146j . which will be different in CS gravity relative to GR. Smith, 
et. al. [37] studied the secular time variation of the longitude of the ascending 
node f2 or b |147j in the non-dynamical modified theory and found it to be given 

f2 or b = ^gr + ^cs, (237) 
where the GR Lense-Thirring drag is given by 

n ° R = aHl-e*y (238) 

with eccentricity e, the magnitude of the spin angular momentum of the central 
body J and the CS correction fics given by 

15a 2 

^cs = -^2~j2(m C sR)yi(m C sa), (239) 

with the semi-major axis a, Earth's radius R, and the spherical Bessel func- 
tions of the first and second kind je(-) and yt{-). 

The LAGEOS satellites have measured Cl and found it in agreement with 
GR up to experimental error, which thus allows for a test of non-dynamical 
CS gravity [37]. Figure [2] shows the ratio of the GR and CS predictions as 
a function of the characteristic CS mass |37J, where the shaded region cor- 
responds to a la deviation from the experimentally measured value. The re- 
gion where the CS prediction is in agreement with experiment is then all 



m cs ^ 2.5 x 10 km i I . A 2a or 3a constraint would increase the shaded 



region by roughly a factor of two or three, thus forcing mcs ^ 2 x 10~ 3 km -1 
or mcs ^ 1.5 x 10~ 3 km~ . Taking the most conservative estimate translates 
into a bound for the CS scalar of < 3000 (n/a) km with more than 99% 
confidence. 

Another type of torque-induced precession is also affected by CS modified 
gravity, namely that experienced by a gyroscope. Consider then a gyroscope 
with spin angular momentum S in circular orbit around the Earth. Neglecting 
geodetic precession, which is unaffected in the non-dynamical modified theory, 
the rate of change S is given by 

gi = 2e oi ^B j S k , (240) 

where B % is the gravitomagnetic field. The precessional angular velocity is then 



10 The relation in Eq. (|239|) is non-monotonic, and thus, many isolated islands exist 
for which the observed precession is consistent with the CS prediction. Nonetheless, 
it is always true that for the range of values quoted here the observed precession is 
always consistent with the CS prediction. The bound quoted by Smith, et. al. are 
somewhat more stringent, thus including more island but also including regions of 
the parameter space that are inconsistent with observations [148 
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Fig. 2. This figure shows the ratio ilcs/^GR as a function of CS characteristic 
mass for the LAGEOS satellites, with semi-major axis a ~ 12, 000 km. The shaded 
region corresponds to a 10% experimental error bar with a la detection confi- 
dence [13011311132] . 

given by $ := | J S , '|/| ( S'*|, which is corrected in CS by 

t — = ^—j 2 (m cs R) [yi(mcsa) + m CS GW/o(m cs a)] , (241) 
$ GR R 2 

where <3>gr is the GR prediction [37]. In Eq. (12411) . R is the distance from the 
center of the Earth to the gyroscope (roughly 7000km for Gravity Probe B) 

Given an experimental verification of the Lense-Thirring gyromagnetic ef- 
fect, one could then place a constraint on the CS scalar, as explained previously 
with the LAGEOS satellites. Gravity Probe B |141|ll42j was designed to mea- 
sured precisely this effect to percent accuracy, but since its launch it has faced 
certain difficulties that might degrade its accuracy |149j . Smith, et. al. studied 
the possibility that Gravity Probe B could place a stronger constraint than 
the LAGEOS satellites, but this was found to not be the case [37]. 

The estimates of Smith, et. al. [37] presented above assume a spherical 
Earth, where its internal structure is neglected. The Earth, however is an 
oblate spheroid with layers of different density. These non-spherical corrections 
seem not to matter since the CS correction is affected by them in the same way 
as GR. Thus, the relative difference between the GR and CS effects remains 
roughly the same. Lastly, the estimates described above assume the Earth is 
not moving in the barycenter frame. In the point particle case, we saw that 
if their velocity is non-vanishing then the CS modification leads to a non- 
boundary correction to the gravitomagnetic field. These corrections remain to 
be studied further. 
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7.2 Binary Pulsar Test 



Non-dynamical (and possible dynamical) CS modified gravity has been 
shown to modify only the gravitomagnetic sector of the metric, which does 
not influence most astrophysical processes. This is particularly true outside the 
Solar system, where stars inside some galaxy will possess randomly oriented 
velocities that will lead to a vanishing averaged CS correction. On cosmological 
scales, the CS effect is also mostly irrelevant since, for example, the equations 
of structure formation are not corrected, except for observations of the cosmic 
microwave background, which we shall discuss in the next chapter. 

Some astrophysical process, however, are CS modified. One example of 
this is the formation of accretion discs around protoplanetary systems, but 
although non-zero the CS correction would be difficult to measure because 
it would be greatly suppressed by the almost negligible compactness of such 
systems. Another, perhaps more interesting astrophysical scenario are double 
binary pulsars, such as PSR J0737 — 3039A/B [150J, which consist of two 
rapidly rotating neutron stars that are in orbit around each other. Neutron 
stars possess a mass of approximately 1.4M , while its radius is on the order 
of 10 km, which implies their compactness is roughly 1/5. Such a large com- 
pactness leads to strong gravitational fields that can be used to test GR to an 
unprecedented level [151J. 

Consider then double binary pulsars and let us model their orbital evolu- 
tion via a geodesic study of a compact object in the background of a rotating, 
homogeneous sphere. Under these assumptions, we can employ the gravito- 
magnetic field found by Smith, et. al. [37] , where the motion is determined by 
the four-acceleration a = —4v x B, v is the velocity of one of the binary com- 
ponents and B is the gra vitomagnetic field of the other [Eqs. (11471) and (11481) ]. 



To leading order in 9 1 I . the CS correction to the gravitomagnetic field is 



B cs = °j- cos[£(r)] [3 - tan £ (d x f) - ($ '■ r) r] , (242) 

with f(r) = 2rn/(9a), c = 15aB/{AKR) sin [£(#)], f = f/r, 3 = J/R 2 and 
(•, x) the Euclidean dot and cross products. 

Let us now parameterize the trajectory with equatorial coordinates, where [152j 

f = [cos u, cos i sin u, sin l sin u] , 
i = [— sin u, cos i cos u, sin i cos u\ , 

n = [0, — sin t, cos 2] , (243) 

are the radial, transverse and normal unit vectors relative to the comoving 
frame in the orbital plane. In Eq. (12431) . 1 is the inclination angle, w is the 

11 Employing 9 as an expansion parameter is not formally valid, since this quantity 
can be dimensional. Corrections are actually proportional to O(0/R) or 0(9 /a), 
both of which can be shown to be much less than unity here. 
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argument of perigee, u = f + w with / the true anomaly and = is the 
right ascension of the ascending node, where here this is co-aligned with the 
x vector [152J. 

The variation of the Keplerian orbital elements is determined by the per- 
turbation equations, which can be obtained by projection of the geodesic accel- 
eration onto the above triad. To leading order, the radial projection a r := a-f 
is the only one CS modified, with [143] 

a^ s = — 4co u d {cos i cos [£ (r)] + sin i cos u sin [£(r )] } . (244) 

The precession of the perigee is given by the perturbation equation w = 



—a r /(n a e) cos/, where the mean motion n = J M/a 3 and e is the eccentric- 
ity. 

We shall consider next only the averaged rate of change of orbital ele- 
ments, and in particular < w >. Averaging over one orbital period |143j 

r T uj r 2 * (l-e 2 ) 3/2 

<oj>: = / -dt= / Co — i '- ~df, (245) 

Jo P Jo 2tt(1 + ecos/) 2 V ; 

during which we assume the pericenter is constant and the motion is Keplerian, 
such that 

u ~ / = n (1 + e cos f) 2 (l — e 2 ^j ^ , 

r = a(l - e 2 ) (1 + ecos f)' 1 . (246) 



This last assumption is justified since in the weak field the CS-corrected motion 
of test particles about any background remains CS unmodified. [99J. 

The CS correction to the averaged rate of change of the perigee is then, 
to linear order in the eccentricity [143 J, 

15 JO. {2kR\ . /2an\ , nlH . 

where X := a sin i is the projected semi- major axis. Note that the limit e — > 
is meaningless because the orbital orientation is ill-defined for circular orbits. 
The scaling found here is consistent with the Solar System results presented in 
Eq. (12391) [37] . The ratio of CS correction to the GR expectation is proportional 
to a 2 0/l£ 3 , since < uj > G r~ J/ a 3 . We see then that a Solar System test will 
naturally be weak since a/R ~ 1, while for binaries a/R ~ O(10 4 ). 

Formally speaking, the sole calculation of perigee precession is not suffi- 
cient to constrain an alternative theory with binary pulsar observations [124j . 
In principle, at least two other quantities must be computed and measured 
in order to break the degeneracy in the determination of the individual com- 
ponent masses. In CS gravity, however, other binary pulsar-relevant quanti- 
ties will not be corrected, because Coo identically vanishes for the canonical 
choice of CS scalar, or the correction becomes subleading, as in the case of 
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the quadrupole gravitational wave emission formula. The only post-Keplerian 
parameters that is CS corrected to leading order is Cj in the non-dynamical 
theory with canonical CS scalar and all parameters can be obtained from [153] . 

Yunes and Spergel [143 J used this calculation to place a strong bound on 
the non-dynamical framework. Using the relevant system parameters of PSR 
J0737 — 3039A/B found in |153j . the CS scalar was constrained to be 

< 4 x 1(T 9 km, -> m cs > 100 meV (248) 

which is 10 11 times stronger than current Solar System constraints. A similar 
analysis in the full, strong-field dynamical formalism is still lacking. In partic- 
ular, the coupling of fermions to the CS term might be relevant for calculations 
involving neutron stars [6]. 

7.3 Galactic Rotation Curves 

The only galactic study that has been currently carried out is related 
to the flat-rotation curves of galaxies [129J. Consider a collection of stars in a 
spiral galaxy and measurements of their orbital velocity v and their distance to 
the galaxy center r, through the shift of spectral lines. The plot of this velocity 
as a function of distance is the so-called galaxy rotation curve, which according 
to Newtonian mechanics should obey a square-root fall-off v oc r -1 / 2 . For a 
large number of such galaxies, with different luminosities, Rubin and Ford 
Jr. |154lll55j have found that the galaxy rotation curve flattens with distance, 
which implies the existence of additional non-visible, or dark, matter. 

Konno, et. al. [129J have attempted to explain the flatness of rotation 
curves through non-dynamical CS modified gravity (a = —Ik and (3 = 0, where 
I is a coupling constant) with a non-canonical CS scalar. We have already 
discussed the form of the metric obtained in the slow-rotation limit, when 
assuming $ oc rcos^ [see eg. Eq. (11691) and Eq. (11711) ]. With this solution, 
Konno, et. al. [129J studied circular, equatorial geodesic motion, and found 
that the orbital angular velocity v := rd<p/dt = —Lf/(rE), where E and 
L and the conserved energy and angular momentum, which for large radius 
becomes 

v = ± ] jf + ^ + 0(Jr, (249) 

and we recall thatC^ is a constant that depends on the spin angular mo- 
mentum. This result is to be contrasted with the Kerr metric, for which 
v = (M/r) 1 / 2 - J/r 2 . 

The solution found by Konno, et. al. is interesting in that it sheds light 
on some of the effects of non-dynamical CS modified gravity on certain ob- 
servables with a non-conventional CS scalar, but before victory can be claimed 
over the rotation curves one must consider the solution more carefully. In do- 
ing so, one discovers that this solution possesses a few drawbacks that render 
it rather unphysical as a true spinning, BH solution. One of main problem is 
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rooted in that the solution was found in the non-dynamical formulation, which 
as we have argued is quite contrived, arbitrary and probably not well-posed. 

One then hopes that an embedding of this solution in the dynamical 
framework can be found, but as we explain below this is impossible. The 
dynamical formulation requires that a well-defined scalar field couple to the 
field equations via a stress-energy tensor. The scalar fields studied in the non- 
dynamical formulation, however, all possess the feature of leading to infinite 
total energy in the scalar field. This issue cannot be bypassed by simply stating 
that the stress-energy contribution is second-order in the slow-rotation param- 
eter, since the energy contribution is not large, but infinite. Thus, the solution 
discussed above is not a self-consistent one in the dynamical framework. 

Moreover, on closer inspection, Yunes and Pretorius [84| have found ad- 
ditional solutions in the non-dynamical formulation for other choices of CS 
scalar [see eg. Eq. (j!73p ]. This family of solutions is better-behaved in that 
they do not contain any logarithmic divergences in the metric, but they lead 
to an orbital angular velocity 



in the far- field limit M/r <C 1, which cannot explain the flat rotation curves. 

At this juncture, one could make the argument that the above analy- 
sis is nothing but evidence that Nature somehow selects the solution found 
by Konno, et. al. , but we shall here argue precisely the opposite. The non- 
dynamical framework does not suggest that either of the two solutions dis- 
cussed here is more valid or less valid than the other. In fact, it is the im- 
mense freedom in the choice of CS scalar the leads to two different observables 
and points at an incompleteness of the framework. This observation, coupled 
with the overconstraining feature of the non-dynamical framework, casts some 
doubt as to validity of suggesting the CS correction as an explanation of the 
flatness in galactic rotation curves. 

7.4 Gravitational Wave Tests 

Several tests have recently been proposed of the modified theory with 
GWs. All such tests have so far concentrated on waves generated by binary 
systems in the early inspiral phase, where CS correction arises due to the 
propagation of the wave (instead of the mechanism). As we reviewed in Sec. |5l 
the main effect of the CS correction on the propagation of GWs is an amplitude 
birefringence, characterized by the parameter ( defined in Eq. f)19ip . Thus, if 
a GW detection can constrain the magnitude of (, one can derive a bound on 
the CS scalar. 

Alexander, et. al. [33] have proposed a GW test of non-dynamical CS 
modified with a generic CS scalar through the space-borne GW detector 
LISA |15b'f 157f 158lll59f 160] . The sources in mind are supermassive black hole 
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binaries located at redshifts z < 30. In principle, in order to determine how 
good of a constraint LISA could place on CS modified gravity, one would 
have to carry out a full covariance (Fisher) matrix analysis, including all har- 
monics in the signal amplitude, since the CS correction affects precisely this 
amplitude. 

One can obtain an order-of-magnitude estimate, however, by making the 
following assumptions. First, let there be a GW detection in two GW detectors, 
such that one can reconstruct both the right- and left-polarized amplitudes. 
Second, let us model the noise as white, with one-sided spectral noise density 
Sq. Third, focus attention on the Fisher matrix related to the parameters 
that affect the amplitude of the GW signal, neglecting the phase parameters, 
namely |161lll62j 




where the observation period is (tj, £/), A 1 = (2), x, C) is a vector of parameters 
that affect the amplitude (the natural logarithm of the luminosity distance to 
the source D = lndz, and the cosine inclination angle x — cost) and m RL is 
the scalar detector response. The ensemble average co-variance is then simply 
Vij = (r -1 ),^., where the diagonal components are a measure of the accuracy 
in the determination of the ii parameter. 

For a GW signal detected on plane (% = 0), the accuracy in the determi- 
nation of C is given by [33J 

_(W ? (252) 

"« - v w-d (252) 

where M = m^/ (mi+mj) 2 is the chirp mass, with binary mass components 
mi^, and p 2 is the squared signal-to-noise ratio 

P 2 = j-J t yt( A l + At), (253) 

with the GW amplitudes. In Eq. (I252p . the noise moments are defined 
via 
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For a binary black hole system with total mass M = mi+m 2 = 10 6 M (1 + z) 1 
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at redshift z, a GW one year before coalescence has a wavelength in the range 
A = (10 2 , 10 4 )cHz _1 . For such a system 



(257) 



lO^xoo /lO~ 40 Hz-^ 1/2 




l + Z-y/l 

m = 3-1 x 10- 40 ( — ^— T )(l+z- VT+~zY , (258 



where h WQ = H o /100 sMpc/km. Such a result implies a la upper bound on £ 
of order 10 -19 for a LISA GW detection at z = 15. With the canonical choice 
of i?, this translates to a constraint on the CS scalar of approximately aft/ k < 
10~ 2 km. Not only is the possible constraint on d five orders of magnitude 
better than with Solar System tests, a GW detection also constrains a different 
sector of the modified theory, since it samples the temporal evolution of the 
CS scalar, instead of its local value. This is because a GW detection really 
constrains the evolution of the CS scalar from the time of emission of the GW 
to its detection on Earth. 

Perhaps a more interesting test of CS modified gravity can be performed 
using GWs emitted by extreme-mass ratio inspirals or binary black hole merg- 
ers [99]. These systems sample the strong-gravitational regime of spacetime, in 
which the CS correction is enhanced, as shown by the CS modified Kerr solu- 
tion [84J . The generation of GWs would then also be CS modified, not only due 
to modifications in the trajectories due to corrections to the background met- 
ric, but also due to the fact that the CS scalar must carry energy-momentum 
away from the system. Even ignoring the latter, Sopuerta and Yunes [99J have 
shown that the background modifications lead to extreme- and intermediate- 
mass ratio inspirals, whose waveforms are sufficiently distinct from their GR 
counterpart to allow for a test of the radiative sector of the dynamical theory 
over a few-month integration period. A more detailed data analysis study is 
currently underway to determine the accuracy to which the theory can be 
tested. 



8 Chern-Simons Cosmology 

After the quintessence model was proposed to account for the acceleration 
seen from Type la supernovae [3], Carroll proposed that the quintessence 
field can generically couple to parity violating terms in the electromagnetic 
sector |163|ll64j . This was used as a constraint on the quintessence dark energy 
models since such parity violating, "birefringent" couplings are generic. Soon 
after, the cosmological study of CS modified gravity was first proposed by 
Lue, Wang and Kamionkowski |134j as a way to search for parity-violating 
effects from the GW sector of the CMB polarization spectrum. This was so 
because, in CS modified gravity, an asymmetry of left and right-handed GWs 
leads to an anomalous cross correlation in the GW power spectrum. Such a 
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pioneering study was then followed up by Saito, et. al. [165J, who improved 
on their calculation. 

The possibility of parity violation in the gravitational sector led Alexan- 
der, Peskin and Sheikh- Jabbari (APS) [2"glE?Tj to propose an explanation for 
the observed matter-antimatter asymmetry in the universe. The lepton asym- 
metry is generated from one simple and generic ingredient: the CS correction 
to the action, which leads to an asymmetric left/right production of GWs, 
and due to the gravitational ABJ anomaly [511152] . chiral leptons are pro- 
duced. APS were able to account for a correct amount of lepton asymmetry 
using the bounds placed on the GW power spectrum amplitude and the scale 
of inflation. In what follows, we shall review the inflationary production of 
GWs with the CS term present. We show how the GW power spectrum is 
modified and apply the GW solutions to the inflationary leptogenesis mecha- 
nism of Alexander, Peskin and Sheikh- Jabbari. Following this, we will review 
the LWK analysis of constraining parity violation in the CMB. 

8. 1 Inflation and the Power Spectrum 

In this chapter we shall consider CS modified gravity in the context of 
inflation. The idea of inflation is quite simple: a period in which the 

scale factor a(t) ~ e ff '' r ''' grows exponentially. During such an epoch, the 
FRW equations govern the dynamics of the scale factor a(t) via 

fa\ 2 8n „ k 2a a 2 k _ . 

- = -^Gp + — and — + — + — = -8vrGp, 259 
\aj 6 a A a cr a z 

where H := - is the Hubble parameter. We see that inflation, a(t) ~ e , 
occurs when there is a fluid of negative pressure p = —p. 

In scalar field theories, this requirement can be obtained if the scalar field 
is slowly rolling down its false vacuum potential, often referred to as slow-roll 
inflation. In what follows, we shall assume slow-roll inflation (see eg. |166] for 
further details) and we shall study the power spectrum of CS birefringent GWs 
during the inflationary epoch. We will then rederive the consistency relation 
between the scalar to tensor ratio and the slow-roll index of inflation following 
Ref. j2B]. 

8.1.1 Gravitational Waves and the Effective Potential 

Consider inflation driven by a pseudoscalar field <j) with a standard kinetic 
term and a potential V(<f>, A). Consider also the Einstein-Hilbert action coupled 
to the gravitational CS term, as given by Eq. ([3]) with the choices a = 2k, 
(3 = and /(</>) a functional of the inflation, which is identified with the CS 
scalar. 

Let us now linearize the action in an FRW background, as explained in 
Sec. [51 The field equations for the GW become essentially Eq. (1177P and, 
concentrating on plane-wave solutions, we find the evolution equation for the 
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phase in a left /right-basis [Eq. (11841) ]. Let us now introduce the quantity z s 
via 

z s ( V , k) = a (rj) Jl^kt (260) 

and the new amplitude ^(rj) defined by /i£ = z s hl_, where s = {R,L}, 
\n t L = ± and the subscript k reminds us of the wavenumber dependance. 
The evolution equation for the phase then becomes a parametric oscillator 
equation for 



z" 



(/4)"+ \k 2 -- K = o- ( 261 ) 



z s 



where primes stand for derivatives with respect to conformal time 77. The effec- 
tive potential z"/z s depends on time, on polarization but also on wavenumber, 
which distinguishes this CS effect from the standard GR case, where the ef- 
fective potential depends on conformal time only. 

Let us study the effective potential in more detail. Using Eqs. (12601) 
and (12611) . this potential is given exactly by 

<_«". m if /a 2 )' X s k (f/a 2 )" (X s k) 2 [(f/a 2 f 2 



z s a l-\°k(f'/a 2 ) 2 l-\ s k(f'/a 2 ) 4 [1 - \ s k (f'/a 2 )] 2 ' 

(262) 

For convenience, we choose specify the CS coupling functional /[</>] via [28|3T] 



f " lQn 2 M 2 pi M P1 ' (263) 

where M P1 = m p J v87r is the reduced Planck mass and N a number that can 
be related to the string scale. In terms of the slow-roll parameters e = —H/H 2 , 
5 = —<f)/(H<f)) and £ = (e — S)/H (an overhead dot means a derivative with 
respect to cosmic time), we have at leading order in these parameters (see also 
Ref. rjS7j) 

a{rj) ~ (-77)" 1_e , <p' ~ -M^VTe. (264) 
From these expressions, one deduces that to leading order 

? = -vmo?^ is? fej ^" (265) 

because % ~ —(1 + e) / 77. 

The equation of motion is still in a rather involved form, but it can be 
simplified further by introducing the characteristic scale kc and G parame- 
ter [281I8T] 

^M%J^<> -Ml)'*- (266 > 
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Let us further introduce the variable x = Qkr]/8 < 0, such that Eq. (12611) 
takes the form 

d 2 // . f64 



+ 



with 
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dx 2 

2 + 3e 

x 2 
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// = 0. 
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x(l-X s x) 4(1-A s a;) 2 ' 

Using slow-roll perturbation theory and the functional relations for / ((f)), the 
equation of motion can be maximally simplified into 



dX 

dr 2 



+ 



1 i0 1 
~4 + 16Y + IT 2 



/£ = 0. 



(269) 



where we have defined r = 16i(l + x)/Q. We recognize Eq. (12691) as the 
Whittaker equation [see eg. Eq. (9.220.1) of Ref. [168J], whose corresponding, 
normalized solution is 



£l e ifcr ?i e _7re / 32 W 



i0/16,O 



16i(l + x) 




(270) 



where W KtfJ ,(z) is the Whittaker function. 

8.1.2 CS Corrected Inflationary Gravitational Waves on Large Scales 

Let us now concentrate on super-horizon scales, ie. x ~ 0, such that we 
can make contact with CMB observables. The effective potentials can then be 
approximated via 



x- 



where the first term gives the standard slow-roll term behavior, while the 
second term represents the birefringent CS correction. With this potential, 
the equation of motion becomes 
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(272) 



Using results from Ref. |169j . we define 



• L 256 

y = l V 1 + ~Q5~ X > K = 



256/0 2 



+ e. 



(273) 



and simplify the equation of motion in the limit 256/0 2 ^> 1, which corre- 
sponds to large scale behavior, such that Eq. f)272p becomes 



d 2 /4 
dy 2 



+ 



1 k (\ 



(274) 
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The exact general solution to this equation is given in terms of Whittaker 
functions 

Mkfa) = Cl{k)W K £ (y) + CI(k)W-^ (-y) , (275) 

where C((k) and Cf (fc) are two constants of integration. 

Initial conditions are determined in the divergence- free region — 1 < x <C 
0, where we assume /x£ has a plane wave behavior |169] 

Cl(k) = -^jf exp , C^k) = , (276) 

where £ F1 is the Planck length and we have used that W Ki (.(y) ~ e~ y ^ 2 y K as 
y — ► +oo. This assumption is equivalent to requiring that non-linear phenom- 
ena occurring near the divergence of the effective potential does not affect 
the initial conditions in the region x > —1. Indeed, such an assumption is 
also made in inflationary cosmology, where the vacuum is assumed to be the 
correct initial state, in spite of the fact that modes of astrophysical interest 
today originate from the trans-Planckian region |17Uj . A possible weakness 
of the above comparison is that, in the case of the trans-Planckian problem 
of inflation |170|171j . one can show that (under certain conditions) the final 
result can be robust to changes in the short distance physics |170|171] . In the 
present context, however, it is more difficult to imagine that the non-linearities 
will not affect the initial conditions. On the other hand, in the absence of a 
second-order calculations and as a first approach to the problem, this seems 
to be quite reasonable. 



8.1.3 Power Spectrum and Tensor-to-Scalar Ratio 

The power spectrum can be calculated either as the two-point correlation 
function in the vacuum state or as a classical spatial average. We take here 
the latter view, since a fully consistent quantum formulation of the present 
theory is not yet available. The power spectrum is then defined as 



(hij {v, x) h* 3 (v, x)) = — / dx hij (rj, x) h lJ (rj, x) 



(277) 



with V — J dx is the total volume. Using the properties of the polarization 
tensor, straightforward calculations show that 
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from which we deduce the power spectrum 

k 3 
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(279) 
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The power spectrum is usually proportional to 2k 3 /it 2 , where here the factor 
of 2 is missing because we have not summed over polarizations. 

The CS corrected power spectrum could be calculated exactly in terms 
of the Whittaker function, but only the large scale behavior is needed so, in 
this regime, one has 
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where T[-} is the Gamma function and ao is the value of the scale factor today. 
Eq. (12801) can be expanded to first order in the slow-roll parameter to obtain 
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where \l/ is related to the derivative of the Gamma function via \I/ = T'/T. 
The amplitude of the CS corrected right-polarization state is reduced while 
the one of the left-polarization state is enhanced. Moreover, at leading order 
in the slow-roll parameter, the spectral index remains unmodified, since = 
dln(/c 3 P^) /dlnfc = — 2e for each polarization state. 

Let us now compute how the tensor to scalar ratio T/S in the modified 
theory. In CS theory, the scalar power spectrum is not modified (see also 
Ref. [172] ) and reads [Ml 
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(283) 



while the tensor power spectrum is given by Eq. (12811) . The T/S ratio is then 
given by 
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where we see that the linear corrections in G has canceled out, and one is left 
with a second-order correction only. Alternatively, we can express the above 
result as the fraction CS correction 
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from which it is clear that the CS correction is not observable, since one has 
assumed here that < 10 -5 (ie. for the divergence of the effective potential 
to be in the trans-Planckian region). 

The super-Hubble power spectrum exhibits two interesting regimes: a 
linear and non-linear one. The non-linear regime occurs when krj ~ O -1 , 
because the effective potential controlling the evolution of the linear pertur- 
bations diverges and linear cosmological perturbation theory becomes invalid. 
This divergence occurs for all modes (ie. for all comoving wavenumber k) but 
at different times. The full non-linear regime has not yet been investigated. 

The linear regime is compatible with the stringy embedding of inflationary 
baryogenesis [27], part of which we discussed in Sec. 13.21 In this context, 
is enhanced, possibly leading to resonant frequencies that could be associated 
with the observed baryon asymmetry. Since is completely determined by the 
string scale and string coupling in a model-independent fashion, one obtains 
a direct link between stringy quantities and CMB anisotropies: 



M 10 is the ten-dimensional fundamental scale and g s is the string coupling. 
For reasonable values of string coupling (ie. weak) and the string scale set to 
10 16 GeV, ~ 10 -2 , but the stringy embedding admits much larger values of 
0, forcing the analysis into the non-linear regime. 

Large values of (eg. > 10~ 5 ) require a non-linear calculation, through 
which one could hope to obtain a significant modification to T/S that might 
lead to an observable CMB signature. Unfortunately, it is precisely in this 
observable regime where technical difficulties have prevented a full analysis of 
the cosmological perturbations. 

8.2 Parity Violation in the CMB 

One of the major issues in the standard model of particle physics is the 
origin of parity violation in the weak interactions. While we know that the 
other gauge interactions respect parity, it may be the case that the there 
is a definite handedness on cosmological scales. The polarization pattern in 
the CMB fluctuations can leave an imprint of parity violation in the early 
universe through a positive measurement of cross correlation functions that 
are not parity invariant. 

The measurement of parity violation from CMB polarization was first dis- 
cussed by Lue, Wang and Kamionkowski |134j . They realized that the presence 
of the CS term naturally leads to a rotation of the plane of polarization as a 




(287) 



where we have used that 
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CMB photon travels to the observer. It was later realized by Alexander [30] 
that gravitational backreaction of parity violating modes can lead to loss of 
power for parity-odd spherical harmonics, which lacking a systematic expla- 
nation, is observed in the CMB for low multipole moments. 

Such considerations can be understood by studying the polarization state 
of light as described through the Stokes parameters. Let us consider a classical 
electromagnetic plane-wave with electric field given by the following compo- 
nents: 

Ei(t) = ai sin(ut — ei) and E 2 (t) = a 2 sm(uit — e 2 ) (289) 

where we assume, for simplicity, that the wave is nearly monochromatic with 
frequency u, such that ai, a 2 , t\, and e 2 only vary on time scales long compared 
to oj~ 1 . The Stokes parameters in the linear polarization basis are then defined 

as 

/ = (( ai ) 2 + (a 2 ) 2 ) , Q ee (( ai ) 2 - (a 2 ) 2 ) , (290) 
U ee (2aia 2 cos 5) , V = (2a 1 a 2 sin 5) , (291) 

where 5 ee e 2 — e\ and the brackets signify a time average over a time long com- 
pared to uj^ 1 . The I parameter measures the intensity of the radiation, while 
the parameters Q, U, and V each carry information about the polarization 
of the radiation. Unpolarized radiation (so-called natural light) is described 
by Q — U — V — 0. The linear polarization of the radiation is encoded in 
Q and U, while the parameter V is a measure of elliptical polarization with 
the special case of circular polarization ocurring when a\ = a 2 and 5 = ±7r/2. 
From here on we will simply refer to V as the measure of circular polarization, 
which is technically correct if Q = 0. 

While I and V are coordinate independent, Q and U depend on the 
orientation of the coordinate system used on the plane orthogonal to the light's 
direction of propagation. Under a rotation of the coordinate system by an angle 
4>, the parameters Q and U transform according to 

Q' = Q cos(20) + U sin(20), and U' = -Q sin(20) + U cos(20), 

while the angle defined by 

* 1 ( U \ 
<P = - arctan — , 

goes to $ — cf) following a rotation by the angle 0. Therefore, Q and U only 
define an orientation of the coordinate system and not a particular direction 
in the plane: after a rotation by tt they are left unchanged. 

Physically, such transformations are simply a manifestation of the oscil- 
latory behavior of the electric field, which indicate that Q and U are part of 
a second-rank symmetric trace-free tensor Pij, i.e. a spin-2 field in the plane 
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orthogonal to the direction of propagation. Such a tensor can be represented 

as 

(P \ 

Pij =\ \, (292) 

[o-p) 

in an ortho normal eigenbasis, where P = (Q 2 + U 2 ) 1 ^ 2 is ususally called the 
magnitude of linear polarization. For example, in two-dimensional, spherical 
polar coordinates (9, 0), the metric is flij = diag(l, sin 2 9) and the polarization 
tensor is 



( Q{n) -U(n) sin6 ^ 



(293) 



^ — U{fi) sin# — Q(h) sin 9j 

where we recall that run over the angular sector only. 

The temperature pattern on the CMB can be expanded in a complete 
orthornmal set of spherical harmonics: 

4 i = l + E E 4m) Y (im)(n) (294) 

J l=l m =-l 

where 

<*<L) = Y Q I dhT(n)Y^ lm) (n) (295) 

are the coefficients of the spherical harmonic decomposition of the tempera- 
ture/polarization map and T is the mean CMB temperature. Likewise, we can 
also expand the polarization tensor in terms of a complete set of orthonormal 
basis functions for symmetric trace-free 2x2 tensors on the 2-sphere, 

3Vn) ' ' 
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E E [^j^mHjW+affcnj^yW]. ( 296 ) 
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where the expansion coefficients are given by 

1 f /.i.,r.;i./.> n 1 



4m) = ¥ J ^(n)Y ( ?4*(n), af lm) = -J dn^{n)Y^\h\ 

(297) 

The basis functions Y^ m ^(h) and Y^v^n) are given in terms of covariant 
derivatives of spherical harmonics by 

Y (fm)ij = N l (*W):*j - g^O'^W**) , (298) 



and 

(lm)ij 2 



Y,? N > 



(Y(lm):ik£ k j + ^(Zm):jfe efe i^ j (299) 

where is the completely antisymmetric tensor on the 2-sphere, a colon in 



an index list stands for covariant differentiation on the 2-sphere, and N 2 = 
2(Z — 2)!/(Z + 2)! is a normalization factor. Since the Y"(z m )'s provide a complete 
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basis for scalar functions on the sphere, the Yn^y and Yn^j tensors provide 
a complete basis for gradient-type (G) and curl-type (C) STF tensors, re- 
spectively. This G/C decomposition is also known as the scalar/pseudo-scalar 
decomposition, which is similar to the tensor spherical harmonic decomposi- 
tion of Sec. [5j 

Integration by parts transforms Eqs. (I297P into integrals over scalar spher- 
ical harmonics and derivatives of the polarization tensor: 

N r 

4m) = yJ dhY$ m) (h) oy^fi), (300) 

AT „ 

4m) = jrj dnY$ m) (h) 3y*(AV, ( 301 ) 

where the second equation uses the fact that e % i ± = 0. Given that T and TV,- are 
real, all of the multipoles must obey the reality condition a&jU = (—l) m an_ m \, 
where X = {T, G, C}. The spherical harmonics Y^ m ) and Yn^j have parity 
(— 1)', but the tensor harmonics Yn^j have parity (— 1)' +1 . 

The two-point statistics of the temperature/polarization map is then 
given via 

Cf X ' = (a? lm) (a? lm) )*), (302) 

where the averaging is over all 21 + 1 values of m and over many realizations 
of the sky. This two-point statistic is thus completely specified by the six 
(TT, GG, CC, TG, TC, and GC) sets of multipole moments. If the temper- 
ature/polarization distribution is parity invariant, then Cj° and Cf c must 
vanish due to the symmetry properties of the G/C tensor spherical harmonics 
under parity transformations. 

Parity conservation, however, is a theoretical bias. Lue, Wang and Kamionkowski [134J 
provided the first time physical scenario where C Z TC = Cp c ^ due to par- 
ity violation in the GW power spectrum of the CMEr^l. This physical sce- 
nario consisted of GWs sourced by the CS interaction term in Eq. ([3]) with 
(a, (3) = (1, 0) and i? = /(</>) some polynomial function of the inflaton field 0. 
As we have discussed, CS modified gravity leads to amplitude birefringence in 
GW propagation, which in turn leads to an excess of left- over right-cicularly 
polarized GWs that lead to a non- vanishing Cf c [134J. 

In order to understand this, consider GWs during the inflationary epoch. 
These waves stretch and become classical at wavelengths on the order of A ~ 

where fi ~ 1/f is some CS energy scale, until they eventually freeze 
as they become comparable to the Hubble radius. When the waves exit this 
radius, the fraction of the accumulated discrepancy between left- and right- 
polarized GWs can be estimated through the index 

e ~ {M p /^){H/Mpf{^/H 2 )\ (303) 

12 Lue et. al noticed that the CS term violates both parity and time reflections. 
Thus, since gravity is insensitive to charge, CPT is conserved. 
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where H is the Hubble scale and /' ~ l//^ 2 - The factor H 2 /4> ~ 10 -5 is 
associated with the amplitude of scalar density perturbations, while H/Mp < 
10~ 6 is related to the amplitude of tensor perturbations [134] . 

Since long wavelength GWs produce temperature anisotropies of curl 
type, an excess of left- over right-polarized GWs produces a nonzero Cj° . 
This is because the multipole coeficients (in m \ will be non-vanishing (right) for 
circularly polarized GW 



t \(5m,2 + 8m-2)Aj(k) even 1 (+), 
a{lm) \-i(5 m , 2 -5 m ^)Af(k) oddl (x), 1 } 

a c = U$m,2 + dm,-2)Af(k) even 1 (x), 
{lm) \-i(5 m , 2 -5 m ^ 2 )Af(k) oddl (+), 1 ) 

where Aj' c are temperature brightness functions and (+, x) stand for plus or 
cross, linear GW polarizations (see eg. |173j ). Likewise the multipole coeffi- 
cients for the gradient component of the CMB polarization are similar, with 
the replacement of A l ' for polarization brightness functions. For a left, cir- 
cularly polarized GW, the sign of the even-/ moments is reversed. The above 
equations allow one to understand why parity is not violated with linearly 
polarized GWs. For example, let us assume that only + modes are present, 
then Cf c by construction. However, a right or left, circularly polarized GW 
possesses both + and x modes, and thus the cross-correction is non- vanishing: 

Cf c = 2(2l + l)- 1 Ar(k)Af(k) (306) 

Lue, Wang and Kamionkowski [134] conclude that the parameter e in 
Eq. (13031) could in principle be measured by a post-Planck experiment with 
a sensitivity of 35fiK, a result that was later confirmed by the more detailed 
study of Saito, et. al. [165]. Such results have aroused interest in the polar- 
ization detection community, pushing them to improve their detection sen- 
sitivities to measure CS-inspired, CMB parity violation. For example, Keat- 
ing et al 1 174] have considered non-vanishing parity violating correlations in- 
duced by a class of telescope-beam systematics, which can mimic the birefrin- 
gence effect. Furthermore, other authors have generalized the parity violation 
cross-correlations to CPT violating cross correlations in the photon sector 
[175111761177117801 . 



8.3 Leptogenesis and the Baryon Asymmetry 



Collider experiments have established a symmetry between matter and 
anti-matter, confirming the prediction of baryon number conservation in the 
Standard Model (SM) of elementary particle physics. In the visible Universe, 
however, there is an excess of matter over antimatter supported by the re- 
cent CMB determinations of the cosmological parameters, in particular by the 
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WMAP experiment pp. In view of this, one of the major puzzles in cosmology 
and particle physics is to understand why and how the matter asymmetry was 
generated during the course of the evolution of the Universe starting from a 
symmetric "soup" of matter and antimatter soon after the Big Bang. In this 
section we will show that if CS gravity is active during the inflationary epoch, 
a novel mechanism of leptogenesis is generic. 

Quantitatively, the baryon asymmetry can be expressed in terms of the 
ratio of baryon density excess to photon density excess [179] 

— = (6.5 ±0.4) x 1(T 10 , (307) 

n 7 

where % = n& — ni is the difference in number density of baryons and an- 
tibaryons and n 7 is the number density of photons. This ratio is time indepen- 
dent, as the evolution of the n& and n 7 with the cosmic Hubble expansion are 
identical. Such a large baryon excess cannot be explained within the SM [180j . 
because baryon number violating interactions are here loop suppressed. The 
only SM source of CP violation in the hadronic sector is in the Dirac phase 
of the CKM mixing matrix, which is not enough to explain the asymmetry of 
Eq. (JSTZD- 

One can map the problem of baryon asymmetry to that of lepton asym- 
metry. This is because the SM weak interactions contain processes, mediated 
by sphalerons (577(2) instantons), which interconvert baryons and leptons and 
are thermally activated at temperatures greater than ITeV. Thus, a baryon 
asymmetry can be generated through the generation of net lepton number at 
high temperature through out-of-equilibrium and CP-asymmetric processes 
[181|I182] . Such scenarios are commonly referred to as leptogenesis. 

Leptogenesis is a valid route to explain the baryon asymmetry, provided 
one can construct a model that fulfills the so-called Sakharov conditions for 
leptons |183j : 

(1) Baryon number violating interactions should be present. 

(2) Charge and parity (CP) should be violated. 

(3) CP and baryon number violating interactions should be active when the 
Universe is out of thermal equilibrium. 

These three requirements constitute model-independent, necessary conditions 
to generate a baryon asymmetry dynamically from symmetric initial condi- 
tions. Within the SM, however, there are no such leptogenesis models, and 
hence one is forced to associate the observed baryon asymmetry to physics 
beyond the SM. 

8. 3. 1 Outline of the mechanism 

A mechanism for the creation of baryon asymmetry, associated with grav- 
itational inflationary fluctuations, was presented by Alexander, Peskin and 
Sheikh- Jabbari (APS) [28]. The key to this mechanism is CS modified gravity, 
where we shall here follow APS and set (a, (3) = (8k, 0), 9 = /(</>) and associate 
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with the inflaton field. Referring back to Sec. I2.4[ the function /(</>) must 
be odd in 0, thus implying that /(</>) is odd under P and CP transformations, 
and allowing the parameterization of Eq. (12631) . 

Let us first spell out how the three Sakharov conditions are realized in 
the APS model of baryon asymmetry, so-called gravi-leptogenesis: 

(1) Lepton number violation is generated here via the triangle anomaly, 
discussed in Sec. 13.11 In the SM, the lepton number current (and hence 
the total fermion number density), has a gravitational anomaly [56J: 

d a J? = ^*RR (308) 
16^^ 

where the lepton number current is given by J/ = J2i=L,R^il a ^i + ^""Vi, 
N = Nl — Nr equals three in the SM, 7° are Dirac gamma matrices, 
and £, v denote lepton and neutrino species respectively. In the SM, the 
anomaly is thus a consequence of an imbalance between left- and right- 
handed leptons. 

(2) CP violation manifests itself naturally in CS modified gravity. Here, 
lepton number is generated due to a non-vanishing vacuum expectation 
value of the Pontryagin density in the evolution equation of the CS scalar, 
which is associated with the inflaton. In turn, this density is generically 
non-vanishing during inflation due to GW perturbations. 

(3) Out-of-equilibrium conditions arise due to the (exponential) growth of 
the background spacetime. Such a growth leads to lepton number pro- 
duction that is naturally out of equilibrium. 

The APS leptogenesis model can be naturally realized if the inflaton field 
is associated with a complex modulus field. In such a case, one must guarantee 
that the inflaton potential is sufficiently flat, such that the slow-roll conditions 
are satisfied. The simplest model of this kind is that of single-field inflation and 
a pseudo-scalar as the inflaton, known as natural inflation. Such a model 
can be expanded to include multiple axions, as in N-inflation models [184J. 
Such models fit into extensions of the SM and in string-inspired inflationary 
models l£ . 

In the remainder of this chapter, we shall use particle physics notation and 
conventions, where h — c — 1 and M P i = (87rG)~ 1/2 ~ 2.44 x 10 18 GeV is the 
reduced Planck mass. In particular, we shall follow [186] and choose a = 8k, 
(3 = 0, and 9 = /(0), although we shall often work with the dimensionless 
functional F(<j>) := f ((f>)/Mp V 

8.3.2 Gravitational Wave Evolution During Inflation 

Although GW solutions have already been discussed in Sec. [5l it is in- 
structive to present the expansion of the Lagrangian to second order in the 
metric perturbation. In the TT gauge, assuming a GW perturbed FRW metric 



For a concise review of string-inspired inflation see [185] 
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and using a left /right basis for the GW polarization, (see Eq. ( I175P ): 



L = -(h L nh R + h R Uh L 
+ ltoF(4>) 





— hu 

dt 2 dtdz 



— h — 

dt R dtdz 



h 92 92 h " 
dtdz L dz 2 L dtdz R , 



+ 0{h 4 



d 2 d 2 N 

hL ~ W lL dmJ lR l 



(309) 



hL ~¥t hL 'dWz hR l 



where t stands for cosmic time dt = a(rj)drj and □ = d 2 + 3Hdt — d 2 /a 2 . As 
is clear from Eq. (13091) . if hi and Hr have the same dispersion relation, *RR 
vanishes, while otherwise "cosmological birefringence" is induced. 

From this Lagrangian, the equations of motion for Jil and Iir become 



Uh, 



-2i—h 



L ) 



Uh 



R 



-2i-k 



R ) 



(310) 



where dots denote time derivatives, and primes denote differentiation of F 
with respect to <fi. The quantity G ~ 2F'H<p/Mp h which with the functional 
form of F{4) becomes 6 = VTe/ '(2n 2 )H 2 / 'M 2 P {K, where e = <p 2 /2 (HM P i)~ 2 is 
the slow-roll parameter of inflation |186j . We have here used the fact that the 
inflaton is purely time- dependent and we have neglected terms proportional 
to by the slow-roll conditions. 

Gravitational birefringence is present in the solution to the equations of 
motion. Let us focus on the positive frequency component of the evolution of 
Hl and adopt a basis in which hi depends on (t, z) only. In terms of conformal 
time, then, the evolution equation for Hl becomes 



d 2 , Id, d 2 , 
— h L - 2-—h L - -—r:h L 

dr] z rj arj dz z 



d 2 

-2iQ——h L 

drjdz 



(311) 



which is a special case of Eq. (11771) . If we ignore for the moment and let 
h L ~ e tkz , this becomes the equation of a spherical Bessel function, for which 
the positive frequency solution is 



h+(k,r]) = e +ik(T,+z) (l~ikrj) . 



(312) 



Let us now peel-off the asymptotic behavior of the solution with non-zero 
B. Let then 

h L = e ikz ■ {-ikrj)e mr, g{rj) (313) 
where ^(77) is a Coulomb wave function that satisfies 



dr) 2 ^ 



2kQ 



0. 



(314) 
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which in turn is the equation of a Schrodinger particle with £ = 1 in a weak 
Coulomb potential. For hi, the Coulomb term is repulsive, while for h,R this 
potential is attractive. Such a fact leads to attenuation of Hl and amplification 
of fiR in the early universe, which is equivalent to the exponential enhance- 
ment/suppression effect discussed in Sec. [5j 

As we shall see, generation of baryon asymmetry is dominated by modes 
at short distances (sub-horizon modes) and at early times. This corresponds to 
the limit krj ^> 1. In this region, we can ignore the potential terms in Eq. (I314|) 
and take the solution to be approximately a plane wave. More explicitly, 

g(r]) = exp[ik(l - Q 2 ) 1/2 r](l + a{rj))] , (315) 



where a(rj) ~ log 77/77. 



8.3.3 The Expectation Value of the Pontryagin Density 

Let us now compute use Eq. (13 13[) to compute the expectation value of 
*RR in an inflationary spacetime. This expectation value is dominated by the 
sub-horizon, quantum part of the GW evolution. Hence, to compute ( *R R) 
we only need the two point (Green's) function (h L hn): 



G(x,t] x ,t') = (h L (x, t)h R (x', t')) 



d 3 k 

(2^)3 



For k parallel to z, the Fourier component G k satisfies Eq. (131 ip with a delta- 
function source 



drf 



1 



d 



2(- + k®)— + k 
77 dn 



5(77-7/). 



(317) 



Let us first consider the case where 9 = 0. The solution to Eq. (13171) is 

N ht(k,v)hR(-k,ri') rj < 7/ 



G k o{r),r)') = < 



(318) 



N h L (k,rj)h^(-k,rj') rf < 77 , 



where K = (H 2 /M Pl )k~ 3 /2, h~l is the complex conjugate of Eq. (13121) . and 
hft, hft are the corresponding solutions of the h R equation. Since here = 0, 
Kft~ = Kl~. The leading-order effect of 6 7^ is to introduce an exponential 
dependence from Eq. (I313p . The generic solution to Eq. (13171) is then 

G k = e- k@T >e +k@r >'G k0 (319) 



for both 77 > 77' and 77 < 77'. 

The Green's function in Eq. (I319p can now be used to contract h L and 
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hji and evaluate the quantum expectation value of *RR. The result is 



(*RR) = - [ 

a J 



d 3 k H 2 



(kr/) 2 ■ k 4 Q 



(320) 



(2tt) 3 2k 3 M 2 



where we have neglected subleading corrections in kr] ^> 1. This expectation 
value is nonzero because inflation is producing a CP asymmetry out of equi- 
librium. The above result and computations seem to depend on the choice of 
vacuum state and the form of the Green's function, but in fact this is not the 
case as one can verify by recalculating Eq. (13201) using a different method, 
such as fermion level-crossing |187j . 

8.3.4 Lepton and Photon number density 

The leptogenesis model can be completed by computing the net lepton/ anti- 
lepton asymmetry generated during inflation. Inserting Eq. (13201) into Eq. f l308f) 
and integrating over the time period of inflation, we obtain 



where ni is the lepton number density. The fc-integral runs over all of mo- 
mentum space, up to the scale ]i = 1/F' at which the effective Lagrangian 
description breaks down, but it is dominated by modes at very short dis- 
tances compared to the horizon scale. The ^-integral is dominated by modes 
at large r] and it represents a compromise between two competing inflationary 
processes: exponential expansion to large k and exponential dilution of the 
generated lepton number through expansion. The dominant contribution to 
lepton generation then arises from modes in the region kr] ^> 1. 
Equation (132 ip can be integrated exactly to return 



The factor (H/M P i) 2 is the magnitude of the GW power spectrum. We should 
stress that the usual GW power spectrum comes from the super-horizon 
modes, while the main contribution to n here comes from sub-horizon modes. 
The factor of G is a measure of the effective CP violation caused by birefringent 
GWs, while H 3 is the inverse horizon size at inflation, giving n appropriate 
units. The factor of (fi/H) 6 is an enhancement over one's first guess due to the 
use of strongly quantum, short distance fluctuations to generate *RR, rather 
than the super-horizon modes which effectively behave classically. 

The significance of the lepton number density can be understood by com- 
paring it to the entropy density of the Universe just after reheating, or equiv- 
alently to the photon number density. Recall that almost all of the entropy of 
the Universe is generated during the reheating time and it is carried by the 
massless degrees of freedom, ie. photons. Let us then employ one of the sim- 
plest (and at the same time most naive) reheating model: instant reheating. 




(321) 




(322) 
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We shall assume that all of the energy of the inflationary phase p = 3H 2 Mp l 
is converted into the heat of a gas of massless particles p = TT 2 g^T A /30 in- 
stantaneously. The entropy of this gas is s = 2ir 2 g )l ,T 3 /45, where T is the 
reheating temperature, g* is the effective number of massless degrees of free- 
dom, s — 1.8 [186j . and we have assumed an adiabatic post-reheating 
evolution. With these assumptions we obtain the photon number density 

n 7 = 1.28g; 3/ \HM Pl f 2 . 

Recalling that the ratio of the present baryon number to the lepton num- 
ber originally generated in leptogenesis is approximately np/ni = 4/11 [181j . 
in this model we obtain 

"* 4.05xlO-W^V /2 em 6 , (323) 



n 7 ' \MpiJ \H 

A less naive approach could be to follow the dilution of and p with the 
expansion of the universe to the end of reheating. The final result is the same 
(see, however, [188j for a comment on this point). With the adiabatic expansion 
assumption, Eq. (1323 j) can be compared directly to the present value of rip/n^ 
given in Eq. (13071) . 

8.3.5 Is Gravi-Leptogenesis a Viable Model? 

In order to answer this question, one must numerically estimate rip/n^. 
Substituting for O, Eq. f!323j) becomes 

Clearly, this ratio depends on five dimensionless parameters: g*, H/Mpi, p/Mpi 
and the slow-roll parameter e and N. 

Some of these parameters are already constrained or given by theoretical 
considerations. For example, within the usual supersymmetric particle physics 
models, g* ~ 100 is a reasonable choice. Moreover, WMAP data, through 
the density perturbation ratio Sp/p (for a single field inflation), leads to an 
upper bound on H/M Pi ratio [lj of H/M Pl < 10~ 4 and e < 0.01, or H < 
10 14 GeV. Finally, the factor X ~ 10 2 — 10 10 is inferred from the string theory 
compactification and is proportional to the square of the four dimensional Mpi 
to the ten dimensional (fundamental) Planck mass [27] (see also the Appendix 
in [27]). If the gravi-leptogenesis model is to be viable, we then have that 

~ 10~ 8 - 10°, (325) 

assuming that H saturates its current bound. 

The physically viable range for the parameter p depends on the details of 
the underlying particle physics model on which our gravi-leptogenesis is based. 
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For example, within the SM plus three heavy right-handed neutrinos (and the 
seesaw mechanism), \x could be of the order of the right-handed neutrino mass, 
in which case fi < 10 12 GeV. If we do not restrict the analysis to the seesaw 
mechanism, /x could be larger, up to perhaps approximately the energy scale 
at which the effective field theory analysis breaks down. Therefore, depending 
on the details of the model 

lO- 6 < JL < 1(T 2 - 10" 1 . (326) 
~ Mpi v ; 

Combining all these considerations, one finds that the gravi-leptogenesis 
model predicts 

^10+5q (JL.) 6 . (327) 
n 7 \Mpi) v ; 

Gravi-leptogenesis thus possesses a large parameter space to explain the ob- 
served baryon asymmetry of the Universe, but its final prediction strongly 
depends on the ultraviolet cutoff of the effective theory |189)190j . In turn, this 
cutoff dependence arises because lepton number is generated by graviton fluc- 
tuations, which we know are non-renormalizable in GR. This issue and that 
of initial conditions are currently being investigated. 



9 Outlook 

We have provided a comprehensive review of CS modified gravity from a 
particle physics, astrophysical, and cosmological perspective. From the particle 
physics point of view, the presence of the CS term leads to a gravitational 
anomaly-cancellation mechanism in the standard model. In cosmology, the 
chiral anomaly works together with inflation to amplify the production of 
leptons, leading to a viable model of leptogenesis. The parity violation present 
in the Pontryagin density leads to direct modifications to GW generation and 
propagation. 

While all of these avenues are promising directions of research, there is 
still much to understand. One important issue that needs to be addressed in 
more detail is that of the coupling strength and the potential of the CS field. 
In string theory, the perturbative evaluation of this coupling suggests a Planck 
suppression, while in LQG it is related to a quantum (Immirizi) ambiguity. 
The potential of the CS field, on the other hand, cannot be evaluated with 
current mathematical methods. Thus, new non-perturbative techniques need 
to be developed to fully evaluate these terms. In lack of a concrete, accurate es- 
timate for these terms, we have taken an agnostic view in this review: to study 
the observational consequences of the CS term and its possible constrains as a 
function of the coupling strength. Such a view allows the possibility to search 
for CS effect, a detection of which would provide theorists with a strong mo- 
tivation to develop the necessary tools to evaluate the CS coupling strength 
and the exact potential of the CS dynamical field. 
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Another issue that remains open is that of backreation and propagation 
effects of fermions in the modified theory. Massive sterile neutrinos may be a 
viable candidate for dark matter. In spiral galaxies, dark matter is expected 
to be at the center of the galaxy where a supermassive black hole might 
reside. In the context of CS gravity, the trajectory of sterile neutrinos would 
be controlled by the modified Dirac equation [6] , which could now be studied in 
the new CS-modified, slow- rotating background of Yunes and Pretorius [M] . A 
modification in the behavior of fermions in highly-dense, strongly-gravitating 
systems could change the rate of black hole formation, such as in the the 
symbiotic supermassive formation mechanism [191J. 

A further avenue of future research that is worth exploring is that of 
an exact solution that represents the exterior gravitational field of rapidly- 
rotating black holes and neutron stars in Dynamical CS modified gravity. 
The slow-rotating solution found by Yunes and Pretorius [M] has provided 
some hints as to what type of modification is introduced to the metric by the 
dynamical theory, but an exact rapidly rotating solution is still missing. Such 
a solution would allow for the study of extreme-mass ratio inspirals and their 
associated GWs in the modified theory, which could lead to interesting GW 
tests via observations with Advanced LIGO or LISA. Moreover, such a solution 
could also be used to study modifications to the location of the innermost 
stable circular orbit and the event horizon, which could have implications in 
X-ray astrophysics. 

A final avenue worth exploring is that of the connection of the dynamical 
formulation and GW theory. The former could serve as a platform upon which 
to test well-defined and well-motivated alternative theories of gravity with 
GW observations. An excellent source for such waves are binary systems, the 
detection of which requires accurate templates. GW tests of dynamical CS 
modified gravity would thus require the construction of such templates during 
both the inspiral and merger. Inspiral templates could be built through the 
post-Newtonian approximation, while the merger phase will probably have to 
be modeled numerically. Such studies would then allow GW observatories to 
explore, for the first time, the non-linear quantum nature of spacetime. 
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